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SECTION  I 


INTRODUCTION 

1 . GENERAL 

In  the  past  decades,  extensive  attempts  have  been  made  to  analyze 
and  predict  response  of  reinforced  concrete  members  and  structural  sys- 
tems. 

A complete  understanding  or  ability  to  predict  this  response  has  yet  to 
be  achieved.  The  heterogeneous  nonlinear  constitution  of  reinforced  con- 
crete challenges  the  analyst  to  predict  micro  behavior. 

The  finite  element  method  provides  the  most  promising  means  for  detailed 
modeling  of  reinforced  concrete.  Some  encouraging  results  have  been  obtained 
using  2-D  finite  element  modeling,  but  the  assumptions  required  are  dissatis- 
fying. 

The  constitutive  behavior  of  plane  concrete  and  the  effect  of  confining 
steel  are  highly  dependent  on  3-dimensional  concrete  material  behavior  and 
steel -concrete  interaction.  The  significance  of  this  study  is  to  attack 
the  problem  with  three-dimensional  nonlinear  modeling. 

A detailed  statement  of  the  specific  purpose  and  scope  of  this  study 
is  given  in  this  section,  along  with  an  extensive  definition  of  terms  and 
a summary  of  the  contents  of  this  volume. 
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2.  PURPOSE  AND  SCOPE 

As  work  progressed  on  the  first  task  of  this  study.  Literature  Review 
(Volume  I),  the  objectives  were  redefined  to  a more  realistic  work  scope 
under  the  existing  time  and  budget  constraints.  The  initial  purpose  and 
scope  are  reviewed  below  followed  by  the  final  definition  of  work  scope 
established  for  this  study. 

The  purpose  of  this  investigation  was  to  more  accurately  define  the 
response  of  reinforced  concrete  structures  under  complex  dynamic  (blast 
type)  loads  to  the  point  of  collapse.  The  intent  was  to  develop  yield  and 
failure  criteria  for  general  structural  shapes  based  on  a consideration 
of  the  three-dimensional  state  of  internal  loading  and  strain.  (Signifi- 
cant terms  such  as  failure  and  yield  are  defined  in  detail  in  subsection 
3.0,  but  briefly,  failure  defines  the  inability  of  a section  to  develop 
or  maintain  appreciable  resistance  to  further  deformation;  yield  is  the 
onset  of  inelastic  behavior.) 

Proposed  alternatives  for  failure  criteria  were  to  be  evaluated  by  numerically 
analyzing  beam-column,  plate,  and  shell  elements  subjected  to  dynamic  loads. 
Comparison  with  actual  experimental  data  was  to  be  used  to  aid  the  selection 
of  the  final  criteria.  Once  the  criteria  were  selected,  they  were  to  be 
used  to  develop  simplified  design  aids  which  in  turn  were  intended  to  be 
used  as  an  improvement  over  the  design  aids  currently  presented  in  the 
Air  Force  Manual  for  Design  and  Analysis  of  Hardened  Structures  (AFWL-TR- 
74-102). 

Through  joint  discussion  and  concurrence  with  AFWL  personnel,  PMB  Systems 
Engineering  (PMB)  and  Professor  A.  C.  Scordelis  (PMB  consultant),  the  final 
work  purpose  and  scope  were  defined  as  follows: 
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The  overall  objective  to  more  accurately  define  the  response  of  rein- 
forced concrete  members  and  structures  under  complex  dynamic  (blast  type) 
loads  to  the  point  of  collapse  remained  as  the  central  purpose.  In  view 
of  the  very  large  loads  associated  with  blast  load  effects,  allowance  for 
inelastic  or  plastic  behavior  is,  of  course,  essential  to  maintain  practical 
structural  proportions.  Additionally,  consideration  of  three-dimensional 
interaction  of  deformation  seems  vital  to  improving  existing  work. 

Thus,  it  was  felt  that  the  most  important  contribution  in  advancing 
the  understanding  of  reinforced  concrete  behavior  when  subjected  to  the 
very  large  deformation  strains  anticipated  in  dynamic  blast  loads  was  to 
employ  a consistent  3-dintensional  analytical  constitutive  model  for 
studying  various  structural  sections  of  interest.  The  effect  of  confining 
steel  on  response  and  the  yield  and  failure  of  plain  concrete  are  highly 
dependent  on  3-diniensional  stress-strain  relationships.  The  objective  then 
was  to  establish  a state-of-the-art  3-D  finite  element  model  accounting  for 
these  important  nonlinear  effects.  Using  this  analytical  tool,  perform  para- 
metric studies  on  two  selected  structural  sections:  (1)  rectangular,  and 
(2)  cylindrical. 

These  parameter  studies  were  then  used  to  evaluate  failure  strength 
of  various  reinforced  concrete  sections  using  a three-dimensional  finite 
element  analysis  (3-D  FEA)  technique  along  with  a 3-D  nonlinear  Constitutive 
Material  Model.  Each  section  was  subjected  to  a fixed  ratio  of  imposed 
displacement  and  rotation  which  was  monotonical ly  increased  till  the  failure 
strength  of  the  section  was  reached.  Each  section  was  subjected  to  seven 
such  deformation  ratios  which  yielded  seven  different  failure  strengths.  These 
failure  strengths  were  then  plotted  in  terms  of  failure  displacement  and  fa'1 'ure 
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rotation  to  give  a section  failure  envelope.  This  failure  envelope  was  then 
normalized  to  dimensionless  form  of  strain  and  curvature. 

Once  a series  of  such  failure  envelopes  were  obtained  for  various  sections, 
they  were  reduced  to  design  failure  envelopes.  For  example,  a rectangular 
section,  having  a longitudinal  reinforcement  of  2 percent  in  the  tension  and 
compression  zones  and  a certain  area  of  transverse  reinforcement,  the  ultimate 
strength  in  uniaxial  compression  of  concrete  was  varied  to  5 ksi,  6 ksi , 9 ksi 
and  9.3  ksi.  Thus  for  this  rectangular  section,  four  design  failure  envelopes 
were  obtained  corresponding  to  the  four  strengths  of  concrete  in  uniaxial 
compression. 

The  failure  strength  of  the  sections  obtained  from  a 3-D  FEA  was  then 
compared  to  the  ultimate  strength  as  obtained  from  one-dimensional  fiber 
analysis  (1-D  FA).  The  increase  in  strength  using  a 3-D  FEA  is  quite  obvious. 
Ductility  of  the  section  using  a 1-D  FA  was  also  compared  to  the  ductility  of 
the  section  using  3D  FEA. 

To  evaluate  the  failure  capacity  of  complete  structural  systems,  analytical 
techniques  will  have  to  be  developed  to  extend  the  failure  strengths  of  the 
sections  to  predict  failure  strength  of  complete  structural  systems.  This  Is 
beyond  the  scope  of  this  report.  However,  a limited  application  technique 
of  the  failure  design  curves  is  addressed  in  Section  IV. 

3.  DEFINITION  OF  TERMS 

It  will  be  most  useful  to  define  the  important  terminology  used  in  the 
study  In  the  context  of  two  different  analysis  techniques:  (1)  fiber 
analysis  (FA)  and  (2)  3-D  Finite  Element  Analysis. 
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a.  Fiber  Analysis  (FA) 

Ultimate  strength  of  a reinforced  concrete  section  can  be  evaluated 
in  terms  of  its  moment  and  axial  force  carrying  capacities  using  a one- 
dimensional analysis  defined  herein  as  the  fiber  analysis  (Figure  1.1). 

In  such  an  approach,  the  section  is  subjected  to  an  imposed  assumed 
strain  distribution  through  its  depth  which  is  resisted  by  internal 
compressive  and  tensile  forces.  Ultimate  strength  is  determined  by 
the  strains  in  the  section. 

The  compressive  forces  are  due  to  the  concrete  block  in  compression, 
and  the  longitudinal  steel  in  compression  if  it  is  a doubly  reinforced 
section.  These  compressive  forces  in  a doubly  reinforced  section  can  be 
evaluated  from  the  stress  distribution  using  uniaxial  stress-strain 
relationships  of  concrete  and  steel.  The  tensile  internal  resisting 
forces  are  due  to  the  straining  of  steel  in  the  tension  ?one.  The  ability 
of  the  concrete  to  resist  tensile  forces  may  or  may  not  be  considered.  In 
case  of  a singly  reinforced  concrete  section  ultimate  strength  of  the 
section  can  be  reached  in  three  different  ways: 

(1)  Tension  Failure 

This  type  of  failure  occurs  when  the  strain  distribution  through 
the  section  is  such  that  the  compressive  strain  in  the  outermost 
concrete  fiber  reaches  ultimate  strain  of  the  uniaxial  stress- 
strain  relationship  of  concrete  (point  7 in  Figure  1.1)  and  the  steel 
tensile  strain  exceeds  yield  (point  8,  Figure  1.1).  Such  a condition 
occurs  in  an  under-reinforced  section. 

(2)  Balanced  Failure 

This  type  of  failure  occurs  when  the  compressive  strain  in  the 
outermost  concrete  fiber  reaches  ultimate  value  (point  7,  Figure  1.1) 
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Compression 


3-D  Finite  Element  Analysis 

Failure  strength  of  a section  Is 
reached  when  every  element  loaded 
In  compression  crushes. 


Fiber  Analysis 

Ultimate  strenoth  of  a section 
Is  reached  when  the  extreme  concrete 
fiber  In  compression  reaches  a value  t 
ultimate  strain. 


Figure  1.1  Definition  of  Terms.  3-D  Finite  Element  Analysis  and  Fiber  Analysis 
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and  the  tensile  steel  strain  reaches  yield  simultaneously  (point  8, 
Figure  1.1).  Such  an  ultimate  strength  condition  is  reached  in  a 
balanced  section. 

(3)  Compression  Failure 

This  type  of  failure  occurs  when  the  tensile  steel  sr  ain  is 
less  than  yield  (i.e.,  the  steel  remains  elastic),  while  the  outermost 
fiber  of  concrete  in  compression  reaches  an  ultimate  strain  value 
(point  7,  Figure  1.1).  Such  condition  occurs  in  over-reinforced 
sections. 

In  all  three  failure  modes  mentioned  above,  ultimate  strength  occurs 
when  the  outermost  concrete  fiber  in  compression  reaches  its  ultimate  strain 
value  obtained  from  uniaxial  compressive  stress-strain  relationship  of 
concrete  (point  7,  Figure  1.1).  This  value  is  usually  assumed  to  be 
0.003  in/in.  The  resulting  internal  resisting  forces  form  a couple  which 
can  then  be  transferred  into  an  equivalent  resisting  moment  and  axial  force 
at  the  centroid  of  the  section.  The  resisting  moment  and  axial  forces 
associated  with  the  strain  distribution  that  causes  the  ultimate  strength 
of  the  section,  is  called  the  ultimate  moment  and  the  ultimate  axial  force 
carrying  capacities.  These  quantities,  along  with  the  corresponding  cen- 
troidal  rotation  and  axial  deformations,  are  classified  as  the  section 
response  quantities. 

The  history  of  these  response  quantities  can  be  traced  from  zero 
strain  to  a strain  distribution  when  ultimate  strength  of  the  section  is 
reached.  The  response  quantities  can  be  graphically  represented  in  terms 
of  the  moment-rotation  and  axial  force  deformation  curves.  Several 
points  on  these  response  curves  can  now  be  defined: 
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(1)  Linear  Behavior 

A section  is  defined  as  linear  when  the  steel  and  the  concrete 
are  elastic.  Steel  in  the  section  is  considered  elastic  if  the 
strain  in  the  reinforcing  steel  is  less  than  its  uniaxial  yield 
strain.  Assuming  that  concrete  does  not  resist  any  tensile  strain, 
the  concrete  section  is  considered  linear  until  the  compressive 
strain  in  the  outermost  fiber  reaches  the  value  of  strain  at 
maximum  stress  of  the  uniaxial  compressive  stress-strain  curve  of 
concrete. 

(2)  Cracking  Point 

This  point  in  the  M-0  and  P-S  response  curves  is  defined  when 
the  strain  distribution  in  the  section  is  such  that  the  concrete 
in  the  tension  zone  has  reached  its  cracking  strength  (point  5, 

Figure  1.1).  This  cracking  strength  may  be  reached  when  the 
stress  in  the  extreme  concrete  fiber  in  tension  reaches  a value 
of  approximately  10%  of  the  ultimate  compressive  strength. 

(3)  Yield  Point 

Yield  point  on  the  M-0,  P-6  response  curves  are  defined  when 
the  strain  distribution  is  such  that  the  tension  zone  steel  reaches 
its  yield  strain  value  or  the  outermost  fiber  concrete  strain  in 
compression  reaches  a value  of  strain  at  maximum  uniaxial  compressive 
stress  (point  6,  Figure  1.1).  Thus,  if  one  has  an  under- rein forced 
section, yield  strength  is  determined  by  tension  zone  steel,  and  for 
balanced  and  over-reinforced  sections,  yield  strengths  are  defined 
by  the  strain  in  the  outermost  fiber  of  concrete  in  compression. 


(3)  Ultimate  Strength 

Ultimate  strength  of  the  section  is  defined  when  the  strain 
in  the  outermost  fiber  of  concrete  in  compression  reaches  an 
ultimate  value  of  uniaxial  strain  (point  7,  Figure  1.1).  Three 
different  common  types  of  ultimate  strengths  have  already  been 
defined  in  preceeding  paragraphs. 

To  suimiarize  the  ultimate  strength  obtained  from  a fiber  analysis: 

1.  The  ultimate  strength  of  the  section  is  reached  when  the 
outermost  fiber  of  concrete  in  compression  reaches  an 
ultimate  value  of  uniaxial  compressive  strain,  usually 
0.003  in/in.  Once  such  a condition  is  reached  the  moment- 
carrying capacity  of  the  section  is  considered  to  be  zero 
on  the  assumption  that  this  failure  propaaates  through 
the  depth  of  the  concrete  section  instantaneously. 

2.  The  position  of  the  neutral  axis  of  the  section  constantly 
changes  through  the  depth  of  the  section  along  the  M-9 
and  P-5  curves  from  zero  to  the  ultimate  strength  of  the 
section. 

3.  Confinement  effects  are  not  considered.  The  member  or 
the  section  is  considered  as  a one-dimensional  element. 

b.  Three-Dimensional  Finite  Element  Analysis 

The  second  approach  to  evaluation  of  the  ultimate  strength  is  to 
use  a three-dimensional  finite  element  analysis  technique  in  conjunction 
with  a 3-D  nonlinear  constitutive  model  (Figure  1.1).  The  ultimate 
strength  will  be  referred  to  as  the  failure  strength  whenever  associated 
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with  a 3-D  finite  element  analysis  as  opposed  to  a fiber  analysis 
evaluation.  This  report  essentially  comprises  of  evaluating  the 
failure  strength  of  various  reinforced  concrete  sections  using  a 
3-D  ideal ization. 

The  first  major  difference  between  the  conventional  fiber  analysis 
technique  and  the  3-D  finite  element  analysis  (3-D  FEA)  technique  adopted 
here  is  the  definition  of  the  ultimate  strength  of  the  fiber  analysis  as 
compared  to  the  failure  strength  of  the  3-D  FEA  technique.  Failure 
strength  in  a 3-D  FEA  is  defined  when  the  equivalent  plastic  strain  in 
every  layer  of  the  concrete  in  the  compression  zone  reaches  a value 
of  equivalent  plastic  strain  for  concrete.  The  equivalent  plastic 
strain  is  evaluated  using  classical  laws  of  plasticity.  The  magnitude 
of  ultimate  plastic  strain  is  evaluated  from  concrete  uniaxial  compressive 
stress-strain  relationships. 

In  a 3-D  FEA,  the  stress  and  strain  is  usually  calculated  at  the 
integration  points  of  the  3-D  finite  elements  as  compared  to  the  fibers 
in  the  one-dimensional  fiber  analysis.  Thus,  if  the  equivalent  plastic 
strain  of  an  integration  point  associated  with  the  outermost  element  of 
concrete  in  compression  is  beyond  the  equivalent  ultimate  plastic  strain, 
that  integration  point  is  assumed  to  have  reached  its  failure  strength. 

It  can  no  longer  resist  any  loads  or  deformations.  The  loads  are  then 
redistributed  to  the  remaining  elements.  If  the  equivalent  plastic 
strains  of  all  integration  points  in  the  compression  zone  are  beyond 
the  equivalent  ultimate  plastic  strain,  then  the  failure  strength  of  the 
section  is  reached. 
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An  integration  point  may  achieve  its  failure  strength  by  three 
different  modes.  These  modes  are  now  defined: 

(1)  Cracking 

When  an  integration  point  is  in  isotropic  tension,  and  the  effec- 
tive stress  at  that  point  as  evaluated  by  the  3-D  constitutive  material 
model  is  equal  to  the  effective  cracking  stress,  the  cracking  strength 
of  the  integration  point  is  said  to  have  been  reached  (point  1, 

Figure  1.1).  The  cracking  strength  of  the  section  is  reached  when 
all  the  integration  points  in  tension  reach  their  cracking  strength. 
This  definition  of  cracking  strength  is  different  from  the  fiber 
analysis,  meaning  where  once  the  cracking  stress  in  the  outermost 
fiber  of  concrete  in  tension  is  reached,  the  section  cracking 
strength  is  defined. 

(2)  Initial  Yield 

Initial  yield  of  the  section  in  3-D  FEA  is  defined  when  the 
effective  stress  at  an  integration  point  exceeds  or  is  equal  to 
the  equivalent  stress  defined  by  the  first  nonlinearity  of  the 
uniaxial  compressive  stress-strain  relationship  (point  2,  Fiqure  1.1). 
Initial  yield  of  the  section  may  be  defined  where  the  outermost 
layer  of  concrete  in  compression  reaches  such  a stress  state.  Note 
the  similarity  in  definition  of  initial  yield  in  the  case  of  3D  FEA 
as  compared  to  1-D  FA. 

(3)  Crushing 

Crushing  of  an  integration  point  in  3-D  finite  element  analysis 
is  defined  in  two  ways,  depending  on  the  stress  state.  In  a tension- 
compression  stress  state,  the  crushing  strength  is  said  to  have  been 
reached  if  the  equivalent  plastic  strain  is  equal  to  or  greater 


than  the  equivalent  plastic  fracture  strain,  corresponding  to  the  strain 
at  ultimate  stress  of  the  uniaxial  compressive  stress-strain  relationship 
of  concrete  and  laws  of  plasticity  (point  3,  Figure  1.1).  If  on  the  other 
hand  the  integration  point  is  in  isotropic  compression,  then  the  crushing 
strength  is  reached  when  the  equivalent  plastic  strain  evaluated  is  greater 
than  the  equivalent  ultimate  plastic  strain  calculated  using  classical 
laws  of  plasticity  corresponding  to  the  ultimate  strain  of  the  uniaxial 
compressive  stress-strain  curve  (point  4,  Figure  1.1). 

(4)  Failure  Strength 

Section  failure  strength  is  assumed  to  have  been  reached  when  every 
integration  point  of  all  the  finite  elements  loaded  in  compression  crush 
either  in  compression-compression  or  tension-compression, 
c.  Comparative  Summary  of  FA  and  3-D  FEA 

The  basic  differences  between  the  one-dimensional  fiber  analysis  and 
the  3-D  FEA  are  as  follows: 

1.  The  ultimate  strength  of  the  reinforced  concrete  section  is 
assumed  to  have  been  reached  in  the  one-dimensional  analysis 
when  the  strain  in  the  outermost  fiber  of  concrete  in  compres- 
sion reaches  the  ultimate  value  of  strain  of  the  uniaxial  com- 
pressive  stress-strain  relation.  In  contrast,  when  using  3-D 
FEA  and  an  associated  3-D  nonlinear  constitutive  material 
model,  failure  strength  of  the  section  is  defined  when  every  3-D 
finite  element  of  concrete  loaded  in  compression  reaches  its 

crushing  strength.  j 

2.  While  evaluating  the  failure  strength  of  a section  by  the 
3-D  FEA  technique,  the  section  was  subjected  to  a mono- 
tonically  increasing  deformation  pattern  in  a fixed  ratio 
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of  displacement  and  rotation.  Thus,  the  neutral  axis  of 
imposed  deformations  was  kept  fixed  while  these  deforma- 
tions were  increased,  until  the  failure  strenqth  of  the 
section  was  reached.  In  contrast,  in  a fiber  analysis 
the  neutral  axis  of  imposed  deformations  was  kept  fixed 
while  these  deformations  were  increased,  until  the  failure 
strength  of  the  section  was  reached.  In  contrast,  in  a 
fiber  analysis  the  neutral  axis  changes  position  within 
the  section  as  defined  by  the  imposed  strain  distributions. 
Hence  in  a fiber  analysis,  all  points  in  the  M-e  response 
curve  have  different  neutral  axis  positions,  while  in  the 
case  of  3-D  FEA,  with  imposed  deformations,  the  neutral 
axis  is  kept  constant  for  all  points  on  the  M-0  and  P^5 
curves,  for  a given  imposed  deformation  ratio. 

3.  No  effect  of  confinement  is  considered  in  the  1-D  analysis 
while  the  3-D  FEA  with  a 3-D  constitutive  material  model, 
does  account  for  the  effect  of  confinement. 

4.  It  is  appropriate  at  this  stage  to  define  ductility  of  a 
section.  The  displacement  or  rotational  ductility  of  the 
section  is  defined  by  the  ratio  of  displacement  or  rotation 
of  the  section  at  ultimate  or  failure  strength  to  displace- 
ment or  rotation  of  the  section  at  initial  yield.  The 
initial  yield  and  ultimate  strength  for  a 1-D  analysis  and 
initial  yield  and  failure  strength  for  a 3-D  analysis  have 
already  been  defined  in  preceding  paragraphs. 
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A detailed  description  of  initial  yield  for  purposes  of 
ductility  evaluation  using  a 3-D  analysis  is  now  presented. 

As  was  mentioned  previously,  the  initial  discontinuous 
surface  or  initial  yield  is  defined  when  the  effective 
stress  evaluated  using  a 3-D  constitutive  model  is  greater 
than  the  effective  stress  corresponding  to  the  point  where 
the  first  nonlinearity  occurs  on  the  concrete  uniaxial 
compressive  o-e  relation.  This  point  is  arbitrary  and 
defined  by  the  user.  To  have  a more  consistent  definition 
of  ductility,  it  was  assumed  that  for  ductility  calculations, 
initial  section  yield  occurs  when  the  effective  stress 
evaluated  using  a 3-D  material  model  is  greater  than  the 
effective  stress  corresponding  to  65%  of  ultimate  uniaxial 
compressive  c-e  relation  calculated  in  a similar  manner. 

4.  SUMMARY  OF  REPORT  CONTENT 

Based  on  the  Volume  1 literature  survey  of  experimental  results  and 
analytical  models,  the  constitutive  relationship  suggested  by  Chen  and 
Chen  (rafs.  1 and  2)  were  adopted  by  PMB  Systems  Engineering  for  the  char- 
act'  rizi.tion  of  various  sections.  Certain  modifications  were  made  to  the 
proposed  model,  which  were  incorporated  in  a general  purpose  nonlinear 
finite  element  analysis  program  called  ANSR  (ref.  3).  Section  II  describes 

1.  Chen,  A.C.T. , and  Chen,  W.F.,  "Constitutive  Relations  for  Concrete," 

Vol.  101,  No.  EM4,  Journal  of  Engineering  Mechanics,  August,  1975. 

2.  Chen,  A.C.T,  and  Chen,  W.F.,  "Constitutive  Equations  and  Punch  Indenta- 
tion of  Concrete,"  Vol.  101,  No.  EM6,  Journal  of  Engineering  Mechanics, 
December,  1975. 

3.  Mondkar,  D.  P. , Powell,  G.  H.,  ANSR-1  General  Purpose  Program  for  Analysis 
of  Nonlinear  Structural  Response,  EERC  Report  No.  75-37,  Earthquake 
Engineering  Research  Center,  University  of  California,  Berkeley,  December, 
1975. 
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these  constitutive  relationships  and  verification  results.  Two  different 


geometrical  sections  were  analyzed:  a rectangular  section  , and  a cylindri- 
cal section.  Results  of  such  analyses  are  presented  in  Section  III.  Appli- 
cations of  the  proposed  model  are  presented  in  Section  IV.  Proposed  failure 
design  curves  derived  from  the  failure  envelopes  of  various  sections  are 
also  presented  in  this  section.  An  example  application  of  the  design  failure 
envelopes  for  a cylindrical  section  is  shown  in  Section  IV.  Assumptions 
and  limitations  of  the  proposed  model  are  also  presented  in  Section  IV. 
Appendix  A contains  a brief  description  of  the  computer  program,  the  user's 
guide,  and  a sample  input  to  the  program  for  a cylindrical  section,  while 
Appendix  B shows  sample  calculations  for  the  ultimate  envelope  using  fiber 
analysis. 


SECTION  II 


ANALYTICAL  MODEL  AND  VERIFICATION  RESULTS 

1 . GENERAL 

To  help  understand  the  mechanical  behavior  of  concrete,  it  is 
possible  to  study  the  response  of  concrete  when  subjected  to  uniaxial 
loading. 

A typical  stress-strain  curve  under  uniaxial  compressive  loading  for 
concrete  is  shown  in  Figure  1.2.  The  curve  is  generally  nonlinear.  However, 
the  curve  is  linear  up  to  stress  states  of  30%  to  60%  of  ultimate  stress. 

If  unloading  occurs  in  this  range,  the  deformations  are  recoverable.  If 
the  load  is  increased  beyond  the  "linear"  range,  permanent  deformations 
occur,  and  the  curve  becomes  nonlinear.  The  nonlinearity  of  the  curve  is 
caused  by  the  development  of  internal  microcracks.  As  the  load  is  increased 
further,  these  cracks  spread  and  interconnect  until  the  ultimate  stress  is 
reached.  Beyond  this  stress  state,  the  curve  descends  with  a negative  modulus 
until  the  concrete  is  completely  crushed. 

The  behavior  of  concrete  when  subjected  to  uniaxial  tension  is  quite 
different  from  the  behavior  under  compressive  loads.  The  stress-strain 
curve  in  tension  is  linear  up  to  the  ultimate  stress,  when  the  concrete 
ruptures  and  cannot  resist  any  further  load.  The  tensile  stress  limit  in 
one  direction  is  affected  by  the  introduction  of  compressive  stresses  in 
other  directions.  For  a constant  tensile  strain  in  one  direction,  the 
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Figure  1.2  Typical  Stress-Strain  Curve  for  Concrete 


27 


2.  CONSTITUTIVE  RELATIONSHIPS 


Chen  and  Chen's  model  (refs.  1 and  2)  considers  concrete  as  an  elastic 
plastic  strain  hardening  and  fracture  material  in  the  isotropic  compression 
region  and  in  the  tension-compression  region.  The  behavior  of  concrete  is 
defined  by  two  relationships,  a loading  function  and  an  incremental  stress- 
strain  law.  The  loading  function  includes  the  initial  discontinuous  surface, 
subsequent  loading  surfaces,  and  the  fracture  surface.  The  initial  discon- 
tinuous surface  can  be  reached  by  elastic  action;  unloading  within  this  sur- 
face causes  no  permanent  deformations.  Straining  beyond  this  surface  causes 
permanent  deformation.  The  fracture  surface  corresponds  to  the  ultimate 
stress  condition.  Ultimate  strain  defines  the  state  of  rupture  where  the 
concrete  can  no  longer  resist  any  load. 

Two  different  but  similar  functions  are  proposed  for  the  loading  surfaces, 
one  for  the  isotropic  compressive  stress  state  and  the  other  for  the  tension- 
compression  stress  state.  Both  these  functions  depend  on  deviatoric  stresses 
and  hydrostatic  pressure  and  are  defined  in  terms  of  I,,  the  first  invariant 
of  stress  state,  and  J„,  the  second  invariant  of  deviatoric  stress. 
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Fig.  2 Initial  Discontinuous,  Subsequent  Loadinq,  and 
Failure  Curves  in  Biaxial  Principal  Stress  Space 
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where  o , o , o are  normal  stresses  in  directions  x,  y,  and  z,  respectively, 
x y z 

and  t , t , t are  the  shear  stresses, 
xy  y z xy 

If  a and  8 in  equations  (1)  and  (2)  above  are  zero,  the  loading  function 
reduces  to  a von  Mises  loading  function. 

Note  that  both  functions  are  similar  except  for  the  sign  of  the  third 
term.  The  constants  a,  8 are  determined  from  the  material  constants  Aq,  t , 
Au,  which  can  be  determined  from  simple  tests.  The  values  of  <2  determine 
the  shape  of  the  loading  surfaces;  this  value  is  found  to  be  3 from  experimen- 
tal results,  which  makes  the  surface  parabolic. 


(3) 


and 


8 = 


A T2 

Vu 


-At2 
u 0 


T2  - T2 

U 0 


(4) 


where  A . t , A , x are  functions  of  the  ultimate  stress  under  uniaxial 
0 0 u u 

compression  f'c  , under  uniaxial  tension  f't  , under  biaxial  compression  f'^  , 

and  the  initial  yield  stresses  under  similar  conditions  f , f.  and  f.  . 

c t be 

Note  that  the  values  of  these  constants  A , t , A , t are  different  in 

o 0 u u 

the  compression  and  the  tension-compression  regions. 

Thus,  for  the  compression  region 


f * - f2 

be  c 


(27, 


bCTc> 


<Tbc  - ’> 


(27' 


be 


1) 


(5) 
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[fc  ' ?bc  <27c  - Tbc» 
^ 3<2fbc  - fc> 


r..  * fj,c  (2  - <U 


u 

v 

c 


be 


3<2fbc  - 1) 


and  for  the  tension-compression  region 


Ao  . <VTt> 


A 

u 

^ ' 


To  2_ 


1 - fi 


f • f. 
c t 


f' 


V_  ft 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


where  ( ) denotes  the  nondimensional ized  quantity  of  the  corresponding 

terms  to  r • 

In  the  loading  function  (equations  1,  2),  if  r = r , the  function 
reduces  to  the  initial  discontinuous  surface  and  if  t = Ty,  the  function 
defines  the  fracture  surface.  Between  these  two  stress  states,  x is  a 
function  of  equivalent  plastic  strain  eP  characterized  by  uniaxial 
compressive  stress-strain  relationship. 
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Thus,  x = r (c*3)  where 


Beyond  the  ultimate  stress  state  x , unloadinq  of  concrete  occurs 
where  the  loading  surface  decreases  while  the  total  equivalent  strain 
increases.  This  state  of  stress  continues  until  the  equivalent  plastic 
strain  reaches  the  equivalent  ultimate  plastic  strain  value.  At  this 
stage  the  concrete  is  assumed  to  be  totally  crushed  and  cannot  resist 
any  load. 

In  the  isotropic  tension  region,  concrete  is  assumed  to  be  a brittle, 
elastic-fracture  material.  In  this  region,  the  initial  discontinuous 
surface  is  assumed  to  be  the  same  as  the  fracture  surface.  If  the 
stress  state  is  located  inside  this  surface,  the  material  is  assumed  to 
behave  elastically;  if  concrete  is  stressed  beyond  this  surface,  it 
is  assumed  to  have  cracked,  and  hence  cannot  resist  any  load.  The  loading 
function  used  is  the  same  as  the  loading  function  of  the  tension-compression 
region. 

Since  the  loading  functions  in  the  compression  region  and  the  tension- 
compression  region  are  different,  it  is  important  to  determine  the  correct 
stress-state  zone,  in  order  to  calculate  the  correct  equivalent  stress.  Two 
different  criteria  to  check  the  stress  state  were  used  by  researchers  at 
Lehigh  University. 

a.  The  first  one  was  used  in  the  Ocean  Thermal  Energy  Conversion  Program 

by  T.  Y.  Chang  and  W.  F.  Chen  (Fritz  Engineering  Lab  Report  No.  414-7)  (ref.  4) 


4.  Chang,  T.  Y.,  and  Chen,  W.  F. , Ocean  Jhermal  Energy  Conversion  Program, 
Report  No.  414.7,  Fritz  Engineering*  laboratory  Report,"  LehTgV  University, 
August,  1976. 


32 


The  compression-compresion  region  was  established  if  the  conditions 
I i < 0 and  (/J2  + ) < 0 were  satisfied,  and  for  the  tension- 

compression  region.  If  I.  < 0 or  (/J,  + — -)  < 0 were  satisfied.  The 

/I 

tension-tension  stress  state  was  established  if  I,  > 0 and 

(/J>  — ) > 0. 

/3  " 

b.  The  other  criterion  was  used  in  program  EDFFEP  for  OTEC  structural 
systems  by  H.  Suzuki  and  W.  F.  Chen  (Fritz  Engineering  Lab  Report 
No.  414-8)(ref.  5).  Here,  the  compression-compression  region  was 
defined  by 

(1)  lx  > T or 

f't 

(2)  -fc  > I,  > fi  and  I.  - m)  > J, 
and  the  tension-compression  region  by 

( 1 ) I j > r and  I j > - f c or 

(2)  -fc  > I > -f'  and  ( •'I  I,  - ^)< 

It  is  believed  that  the  stress  state  region  is  independent  of  the  stress 
parameters  related  to  any  given  material.  Hence,  the  second  criterion  is  not 
recommended.  A discussion  of  the  first  criterion  is  presented  later  in  this 
Section. 

Once  the  loading  function  has  been  defined,  the  incremental  stress-strain 
relationships  can  be  derived  directly  from  the  loadinq  surface  using  the  normality 
condition.  In  the  three-dimensional  case,  the  incremental  relationship  is 

5.  Suzuki,  H. , and  Chen,  W.  F.,  Ocean  Thermal  Energy  Conversion  Program, 


n.,  ana  cnen,  w.  r.,  ucean  inermai  Energy  conversion  Program, 
Report  No.  414.8,  Fritz  Engineering  Laboratory,  Lehigh  University, 
August,  1976. 
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where 


3 - <2 

n = jg- — in  the  isotropic  compression  reqion 

-3  - <2 

n = — jg in  the  tension-compression  region 

All  parameters  can  be  easily  calculated,  excepting  the  value  of  H. 

The  value  of  H,  defined  as  the  strain  hardening  rate,  can  be  written  in  terms 
of  the  slope  of  the  equivalent  stress-strain  curve: 

H = 2oe  • H1  (16) 

where 

H'  = do  /deP 
e 

a = the  equivalent  stress  = / f(o. . ) defined  in  equations  (1)  and 

' J 

(2) 

ep  = equivalent  plastic  strain  = dep  = (dePj  * 

The  og  ~ cp  curves  in  the  compression  region  and  the  tension-compression 
region  are  different.  These  nonlinear  curves  are  assumed  to  be  identical 
in  form  and  are  obtained  from  tests. 
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Since  the  derivation  of  the  incremental  stress-strain  equation  was 
based  on  the  elastic-plastic  model,  any  finite  element  program  for  elastic- 
plastic  analysis  may  incorporate  this  material  model  without  much  difficulty. 

As  mentioned  before,  the  initial  discontinuous  surface  can  be  reached 
by  elastic  action  only.  To  check  if  this  initial  discontinuous  surface 
had  been  reached  during  loading,  certain  yield  criteria  are  checked. 

At  each  increment  of  strain,  the  elastic  incremental  stresses  are  calcu- 
lated using  the  Incremental  stress-strain  law.  The  yield  criterion  for  the 
element  is  then  checked  at  this  stage.  Thus,  if  f(o^)  - t2  < 0,  the 
element  is  under  elastic  loading,  and  the  stress  increments  are  calculated 
from  the  elastic  incremental  stress-strain  law.  If,  however,  f(a.j)  - t2  > 0, 
then  the  element  is  under  plastic  loading.  In  this  case,  then,  the  elastic 
portion  of  the  stress  increment  is  scaled  down  by  a factor  R which  is 


calculated  as  follows: 


f<5ij  + Rd"ij5  - rJ  ’ 0 


where 


a..  = Total  stresses  at  previous  load  step 
■ 0 

da..  = Stress  increment  at  current  load  increment  calculated 

* vl 

from  elastic  stress-strain  law. 


This  condition  can  hp  satisfied  if 


AR2  + 2BR  + C = 0 


where 
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1 2 ) + no2  + (82  + ax2)  - x2  (21) 


Thus  R = ~B-t-  C-  " A- 


(22) 


Note  n = 


3 - k2  . , 


18 


in  isotropic  compression  region 


-3  - *c2  . 

n = — |g in  tension  compression  region 


The  minimum  positive  value  is  the  required  scaling  factor. 

The  remaining  portion  of  the  incremental  strain  is  given  by 


O - R>  dcfj 


Once  the  elastic  portion  of  the  strain  and  the  inelastic  portion  of 
the  strain  increment  were  determined,  the  incremental  stresses  were 
reevaluated  by  forming  the  new  constitutive  matrix  of  equation  (14).  The 
equivalent  stresses  were  then  evaluated  using  appropriate  loading  functions 
after  stresses  and  strains  were  updated. 

During  the  preliminary  analyses,  it  was  found  that  the  main  difficulty 
of  the  model  was  the  criteria  used  to  establish  the  various  stress  states. 
The  criteria  used  to  check  these  states  were  the  first  invariant  of  stress 
state  Ij,  and  the  second  Invariant  of  the  deviatorlc  stress  J2. 

Thus,  if 


I.  < 0 and  /J  + -i-  < 0 

1 2 ^ 
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then  the  stress  state  that  existed  was  in  the  isotropic  compressive  zone. 
These  criteria  were  found  to  be  inconsistent  with  the  loading  function, 
can  also  be  shown  as  follows: 


= T 


1 - a/3  I, 

The  equivalent  stress  is  then  defined  as 


(23) 


a-  ■ /T!°:uT 


(24) 


For  ag  to  be  a real  value,  the  function  ffc^.j)  must  be  positive. 

However,  in  certain  cases  it  was  found  that  f(a..)  was  negative.  This 

1 J 

negative  value  could  be  possible  if  either  (1  - ^ It)  was  negative  while 
the  numerator  of  the  expression  for  f(a. .)  was  positive  or  vice  versa. 

* vJ 

In  most  cases  the  value  of  this  numerator  is  positive.  Hence,  f(a..) 

* J 

is  usually  positive  if  1 - ^ Ix  is  positive,  and  f(a^)  is  usually 
negative  if  the  value  of  (1  - ^ Ij  is  negative.  Let  us  now  examine 
the  condition  under  which  this  would  be  a negative  value. 

In  the  tension  zone: 


A - A 

n _ u o 

ct  — 2 j- 

T _ T 
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where 


Au  _ fc  ■ f,t 
fc  ' 2 


(26) 
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fc  " 2 


(27) 
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Thi 








,V  . T'c  T't 
(^>  6 

(28) 

(— )J  - ^ 

(fc  6 

(29) 

Substituting  the  values  of  A A r T <n  „ 

V o’  Tu’  o in  the  expression  for  a in 

equation  25,  we  get  a as  a function  of  f^,  f'  , f^t  f 


Thus 


a = 


3 (7'c  - - <?r  - M 

C (f'c  Tt  - fc  ft) 


Thus,  for  1 - | I,  to  be  greeter  then  zero,  we  must  have  1 > 2 I,  , or 

I - k (f'c  ' f't>  - (f'r  ' V 


C (fc  f't  * Vt 


< 1 


If  we  now  have  a stress  state  where 

°n  = f't  . o22  = o33  = t12  = t23  = t31  = 0 


then 


or 


£ (T'c  - V - <7r  - V 
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ITc  T't  - W 


< 1 
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This  condition  can  be  achieved  since  in  most  cases  f't  is  about  10 
percent  of  f'  so  that  the  value  of  the  expression  in  the  brackets  can  be 
guessed  to  be  approximately  10.  Now,  if  we  somehow  have  a stress  state 
where  Ij/f'c  is  greater  than  f't/f'c,  it  may  result  in  a negative  value  of 
f(c. .)  so  that  the  effective  stress  is  a complex  number.  Such  a value  of 

■ J 

effective  stress  has  no  physical  meaning.  However,  it  has  been  found  that  if 
the  effective  stress  (assuming  a tension-compression  region)  is  evaluated 
as  a negative  quantity,  and  we  use  the  material  constants  of  a and  3 of  the 
compression  region  and  then  reevaluate  the  effective  stress,  this  value 
turns  out  to  be  a positive  value,  hence  the  effective  stress  is  real.  This 
could  imply  that  the  criteria  used  to  check  the  various  stress  states  may 
not  be  sufficient  and  may  be  inconsistent  with  the  loading  function. 

Hence,  a new  modification  was  introduced  in  the  program:  if  the  stress 
state  was  established  by  the  condition  of  Ij  , and  (Ij  and  Jj  was  satisfied 
or  otherwise,  and  the  effective  stress  so  calculated  was  found  to  be  a com- 
plex number,  then  the  stress  state  was  switched  to  the  opposite  one.  This 
was  done  for  the  compression  region  and  the  compression-tension  region 
only.  Thus,  if  the  stress  state  established  by  the  criteria  proposed  by 
Chen  and  Chen  indicated  the  stress  to  be  in  the  compression-tension  region 
and  the  effective  stress  calculated  using  the  loading  function  of  that 
region  was  found  to  be  negative,  then  it  was  assumed  that  the  stress  state 
was  actually  in  the  compression  region.  Thus,  the  material  constants  for 
this  region  were  reevaluated  and  the  effective  stress  calculated  using  the 
appropriate  loading  function.  A similar  course  of  logic  would  be  applied 
to  the  compression  region,  and  it  can  be  shown  that  the  loading  function 
is  negative  for  certain  stress  states  which  could  be  avoided  if  the  stress 
states  were  switched  as  mentioned  above. 


To  summarize,  there  are  three  different  stress  states  possible: 
a.  Compression-compression  stress  state  is  established  if  I < 0 and 

/J2  + < 0.  Once  this  state  has  been  established,  the  program  sets 

the  parameter  identified  as  IDTC  as  1.  The  yield  code  printed  with  the 
results  of  the  element  is  the  product  of  IDTC  and  the  parameter  defined 
KODY.  The  parameter  KODY  defines  the  nonlinear  behavior  of  concrete 
in  this  region  as  shown  in  Figure  3,  a plot  of  equivalent  stress  and 
strain.  Thus,  if  the  equivalent  stress  og  < tq,  where  tq  is  defined  in 
equation  7,  KODY  = 0,  thereby  establishing  that  the  stress  state  is 
elastic.  If  e<  e^,  and  the  equivalent  stress  og  is  less  than  tu 
(as  defined  in  equation  8),  then  the  behavior  is  plastic  and  KODY  is 
then  set  to  1.  If  the  equivalent  stress  is  greater  than  tu  and  the 
equivalent  plastic  strain  is  greater  than  or  equal  to  equivalent  plastic 
fracture  strain,  then  the  fracture  state  is  assumed  to  have  been  reached. 
KODY  is  now  set  to  a value  of  2.  If  the  equivalent  plastic  strain  is 
greater  than  equivalent  plastic  ultimate  strain,  the  concrete  is  assumed 
to  have  crushed  and  KODY  is  set  to  3. 


b.  Tension-compression  region  is  established  if  Ij  > 


0 or  /T  + 41:  >0. 


c. 


2 /3 

Once  this  state  is  established,  the  program  sets  the  parameter  IDTC 

to  -1.  The  nonlinear  behavior  of  concrete  in  this  region  is  shown  in 

Figure  4.  Note  that  for  this  region,  the  descending  portion  of  the 

stress-strain  curve  is  absent.  The  values  of  KODY  = 0 and  1 have  the 

same  significance  as  for  the  isotropic  compression  region.  For  this 

region,  however,  KODY  = 2 represents  the  crushing  of  concrete. 

_ I 

Tension-tension  region  is  established  if  It  _>  0 and  /J  + > 0. 

In  this  case,  the  parameter  IDTC  is  set  to  -6.  The  stress-strain 
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mm 


cf  CU 


If  ae  < To  - elastic  state  KODY  = 0 

(c  < c^):  Tu  > oe  l To  " elastic  plastic  state  KODY  = 1 

oe  > xu  - fracture  state  KODY  = 2 
e > e u - crushing  state  KODY  = 3 

Note:  x0  and  tu  values  evaluated  from  equations  7 and  8 for  the 
compression-compression  region. 

d£  = (dCg-dEy)  h 

Fig.  3 Criterion  for  Different  Stress  States 


in  Isotropic  Compression  Region 


I I 

e0  Ef 

If  °e  < To  - elastic  state  KODY  = 0 

e < ef:  tu  > oe  > x0  - elastic  state  KODY  = 1 

oe  > tu  or  c > ef  - crushing  state  KODY  = 2 

N°te:  tq  and  tu  are  evaluated  from  equations  11  and  12  for 
tension-compression  region. 


Fig.  4 Criterion  for  Different  Stress  States 
in  Tension-Compression  Region 
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KODY  = 0 
KODY  = 1 


If  ae  < t0  elastic  state 
cre  > To  cracking  state 

Note:  xQ  used  is  the  same  as  for  tension-compression  region 


Fig.  5 Criterion  for  Different  Stress  States 
in  Tension-Tension  Region 
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curve  for  this  region  is  indicated  in  Figure  5.  As  can  be  seen  from 
this  figure,  there  are  only  two  states  possible;  the  elastic  state 
where  KODY  = 0,  and  the  cracking  state  where  KODY  = 1.  The  cracking 
state  occurs  once  the  equivalent  stress  is  greater  than  or  equal  to 
tq  [as  defined  by  equation  (11)]. 

The  material  model,  together  with  a three-dimensional  8 to  20  node 
finite  element  subroutine,  was  added  to  program  ANSR  (ref.  3),  a general 
purpose  three-dimensional  nonlinear  analysis  program.  The  material  model 
was  essentially  incorporated  in  three  basic  subroutines: 

a.  SUBROUTINE  MATP45,  which  essentially  read  in  the  necessary  material 

data,  evaluated  the  material  constants  Aq,  A^,  t , for  the  compression 
and  tension  regions,  and  also  evaluated  the  constants  a and  £ using  the 
above  material  constants.  Equivalent  plastic  fraction  strain  and 
equivalent  plastic  ultimate  strain  are  also  evaluated  in  this  subroutine. 

b.  SUBROUTINE  CMAT45,  which  evaluated  the  elasticity  matrix  defining 

the  incremental  stress-strain  relationship. 

c.  SUBROUTINE  MRES45,  which  evaluated  the  response. 

These  three  basic  subroutines  in  conjunction  with  a few  other  minor 
routines  were  used  to  define  the  material  behavior  of  concrete  using  the 
constitutive  relationships  mentioned  above.  The  main  program  ANSR  was  also 
modified  to  impose  displacements  and  rotation  on  a section  instead  of  imposing 
forces.  The  imposed  displacements  could  be  applied  in  either  of  the  direc- 
tions of  the  coordinate  axes  x,  y,  and  z or  in  all  three  directions,  or 
any  combination  thereof. 


3.  NUMERICAL  STABILITY  OF  SOLUTION  PROCESS 

The  program  ANSR  to  which  the  elastic-plastic  failure  strain  hardening 
material  model  for  concrete  was  incorporated,  has  various  options  of  numerical 
solution  schemes.  The  two  basic  schemes  available  are  the  Newton-Raphson 
scheme  and  the  constant  stiffness  iteration  scheme.  The  solution  scheme 
utilized  was  a mix  of  the  two  techniques.  Within  each  load  step,  iteration 
cycles  were  prescribed  to  allow  for  stiffness  reformation,  and  constant  stiff- 
ness iterations  were  employed  within  in  each  cycle  to  ensure  convergence  of 
the  solution. 

Concrete  is  modeled  as  a mesh  of  discrete  three-dimensional  8-to-20  node 
finite  elements.  To  ensure  the  stability  of  the  solution  scheme,  several 
steps  had  to  be  taken.  Suppose  a reinforced  concrete  section  is  modeled  as 
shown  in  Figure  6.1  using  concrete  solid  elements  and  three-dimensional  truss 
elements  for  reinforcing  steel.  If  this  section  is  subjected  to  an  imposed 
deformation  pattern  of  extension  and  rotation  in  a certain  ratio,  and  this 
constant  deformation  ratio  of  6/0  is  increased  in  a number  of  load  steps, 
a stress  state  will  be  created  such  that  the  bottom  four  gauss  points  of 
the  bottom  element  in  the  tensile  zone  crack.  Once  this  stress  state  occurs, 
the  stiffness  contribution  of  these  four  gauss  points  is  reduced  to  a percentage 
of  the  original  stiffness  (say  10%)  for  the  next  stiffness  reformation.  If 
in  the  next  load  step  the  top  four  gauss  points  of  this  same  element  also 
crack,  then  the  stiffness  of  this  element  is  set  to  zero.  At  selected  degrees 
of  freedom  the  steel  truss  bar  elements  do  not  contribute  any  stiffness,  so 
that  the  resulting  stiffness  would  be  zero.  The  system  of  equations 
K r = R^  is  solved  using  Gaussian  elimination.  However,  during  the  fac- 
torization of  K into  l D LT,  the  elements  of  the  D diagonal  are  checked  for 
small  numbers  resulting  from  the  zero  stiffness  of  cracked  elements.  Such  a 
(a)  Bar  below  letters  indicates  matrices. 
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Transverse  Steel 


Gauss  Integration  Points 


3-D  Brick  Elements  To  Model  Concrete 


Note:  Concrete  Cover  Was  Not  Considered 


Figure  6.1  Reinforced  Concrete  Section.  Finite  Element  Model 


condition  could  cause  large  nodal  displacements  resulting  in  large  displace- 
ments and  imminent  numerical  instability.  Thus,  if  during  the  factorization 
of  K the  element  of  the  D matrix  produces  numbers  smaller  than  10~s,  the 


diagonal  element  was  replaced  with  a very  large  number  of  the  order  of  10 


1 V 


to  stabilize  the  solution.  Note  that  both  10  S and  1014  are  arbitrary  num- 


bers. Thus,  to  summarize,  the  following  steps  ensure  numerical  stability: 

a.  If  a qauss  point  in  an  element  cracked  or  crushed,  then  the  stiffness 
contribution  of  this  gauss  point  to  the  element  stiffness  is  reduced  to 

a certain  percentage  of  its  original  stiffness.  This  number  is  specified 
by  the  user  as  an  input  to  the  program.  A value  of  10  percent  is 
recommended. 

b.  After  the  element  stiffnesses  are  assembled  into  the  structure 
stiffness  matrix,  the  diagonals  produced  during  the  factorization 

are  checked  for  small  numbers.  If  any  element  in  the  diagonal  is  found 
to  be  less  than  10~s,  this  element  is  replaced  by  a large  number  of  the 
order  1014. 

4.  VERIFICATION  RESULTS 

The  material  model  incorporated  into  program  ANSR  was  checked  by 
comparing  the  analytical  results  against  experimental  results.  Since  Chen 
and  Chen  (refs.  1 and  2)  had  already  checked  their  constitutive  relationships 
and  pr  o;  ted  their  results  compared  to  the  experimental  results  of  Kupfer, 
et  al.,  ^*ef.  6)  it  was  decided  to  check  the  model  against  the  results  of 
Chen  and  Chen  and  Kupfer,  et  al. 

hence,  a series  of  runs  was  carried  out  using  a single  brick  element 
of  concrete  7.9  ins.  by  7.9  ins.  by  2 ins.  (20  cm  by  20  cm  by  5 cm).  The 


6.  Kupfer,  H. , Hilsdorf,  H.  K. , and  Rush,  H. , "Behavior  of  Concrete  under 
Biaxial  Stresses,"  Vol.  66,  No.  8,  pp.  656-666,  Proceedings,  Journal 
of  American  Concrete  Institute,  August,  1969. 
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uniaxial  compressive  stress-strain  curve  used  is  shown  In  Figure  6.2. 

This  brick  was  subjected  to  biaxial  forces  in  a ratio  of  Oj/a,,  the  principal 
stresses  in  two  directions.  This  ratio  was  varied  for  biaxial  compressive 
forces,  i.e.,  ratios  of  o,/a2  of  -1/-1,  -1/0,  -1/-0.522.  -1/-0.0001,  and 
for  tension  compression  forces  of  -1/0.103,  -1/0.204,  and  -1/0.052.  The 
results  obtained  from  these  analyses  were  plotted  against  the  experimental 
results  of  Kupfer,  et  al  and  the  analytical  results  of  Chen  and  Chen.  The 
results  are  presented  in  Figures  7 through  14. 

Figure  7,  which  plots  o,  versus  for  Oj /o2  of  -1/-1,  shows  that  the 
present  model  has  more  stiffness  than  the  Chen  and  Chen  or  Kupfer,  et  al. 
results.  In  this  analysis  the  uniaxial  stress-strain  curve  input  to  the 
program  was  assumed  to  be  linear  up  to  a stress  of  30  percent  of  the 
ultimate  stress  f 1 . Chen  and  Chen  (ref.  2)  have  indicated  in  their  report 
that  the  normality  condition  was  best  satisfied  using  a linear  uniaxial 
stress-strain  relationship  up  to  a stress  euqal  to  60  percent  of  the 
ultimate  stress.  Thus  this  relevant  parameter  was  changed  and  the  results 
of  the  new  analysis  are  indicated  in  Figure  8.  It  can  still  be  observed 
that  the  present  model  response  has  more  stiffness  associated  with  it  than 
the  results  of  Kupfer,  et  al.  and  Chen  and  Chen. 

Figure  9 indicates  the  results  of  the  applied  forces  a, /a,  in  the  ratio 
of  -1/-0.522.  The  correlation  in  this  case  is  much  better  than  for  equal 
biaxial  imposed  forces.  The  correlation  in  this  case  is  much  better  than 
for  equal  biaxial  imposed  forces.  Figures  10  and  11  Indicate  results  for 
Oj/o,  of  -1/-0.0001  where  the  change  in  the  modeling  between  the  two 
analyses  is  the  value  of  initial  yield  point  of  uniaxial  stress-strain 
relationships  of  concrete. 
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Fig.  10  Analytical  Results.  Biaxial  Compressive  Forces  Imposed. 
a,/o,  - -1/-0.0001.  fc/f’c  - 0.60 


The  correlation  in  the  case  of  tension-compression  applied  forces 
seems  to  be  better  than  in  the  biaxial  compressive  applied  forces. 

Figures  12  through  14  indicate  the  results  of  the  analysis  for  various 
ratios  of  loading. 

The  general  conclusion  that  can  be  drawn  from  the  above  analysis  is 
that  the  present  model  has  more  stiffness  associated  with  it.  This  stiff- 
ness is  sensitive  to  the  type  of  uniaxial  stress-strain  relationship  input 
to  the  model.  The  exact  parameters  for  the  uniaxial  a-e  curve  used 
by  Chen  and  Chen  in  their  analysis  were  not  known.  In  this  case,  the 
curve  assumed  may  be  different,  which  can  account  for  the  difference  in 
stiffness. 
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SECTION  III 
PARAMETRIC  DATA 

1 . GENERAL 

Once  the  model  was  tested  and  the  results  indicated  in  the  preliminary 
analysis  found  to  be  adequate,  the  analysis  of  two  types  of  sections  was 
considered:  the  rectangular  section  and  the  cylindrical  section.  The 
analysis  consisted  of  imposing  monotonically  increasing  values  of  a specific 
ratio  of  axial  to  rotational  deformation  to  a given  section.  After  such  an 
increment  of  deformation,  total  stresses  in  the  concrete  and  steel  were  com- 
puted to  determine  if  the  concrete  had  cracked,  yielded,  crushed  or  failed, 
or  the  steel  had  yielded. 

Stresses  for  the  concrete  were  calculated  at  the  integration  points 
of  the  three-dimensional  finite  elements.  These  stresses  were  integrated 
over  the  volume  of  the  section  to  give  equivalent  loads  at  the  nodes  of  the 
element.  The  generalized  axial  force  and  bending  moment  on  the  section  were 
also  calculated  using  these  nodal  loads  and  the  loads  in  the  steel. 

For  concrete,  yield  of  an  integration  point  was  defined  when  the  effec- 
tive stress  at  that  point  (equation  1 and  2)  exceeded  the  effective  yield 
stress.  The  effective  yield  stress  t0  was  calculated  using  equation  7 for 
an  isotropic  compression  state  of  stress  and  by  equation  11  for  a tension- 
compression  state  of  stress. 

For  an  isotropic  compression  state  of  stress  at  an  integration  point, 
crushing  was  defined  when  the  effective  stress  (equation  1)  exceeded  the 
effective  ultimate  stress  tu  (equation  8)  and  the  equivalent  plastic  strain 
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at  that  point  exceeded  the  equivalent  ultimate  plastic  strain.  The  equiva- 
lent plastic  strain  is  defined  by  equation  13  while  the  equivalent  ultimate 
plastic  strain  is  calculated  using  laws  of  plasticity  corresponding  to  the 
ultimate  strain  of  the  uniaxial  compressive  stress-strain  curve. 

For  the  tension-compression  state  of  stress,  crushing  of  the  integration 
point  occurs  when  the  effective  stress  of  that  point  (equation  2)  exceeds  the 
effective  ultimate  stress  (equation  12)  and  the  equivalent  plastic  strain 
(equation  13)  exceeds  the  equivalent  fracture  plastic  strain.  The  equivalent 
fracture  plastic  strain  is  calculated  using  laws  of  plasticity  corresponding 
to  the  strain  at  ultimate  stress  of  the  uniaxial  compression  stress-strain 
relation. 

Cracking  of  an  integration  point  in  isotropic  tension  is  defined  when 
the  effective  stress  of  that  point  (equation  2)  exceeds  the  effective  cracking 
stress  (equation  11). 

Once  an  integration  point  of  the  three-dimensional  finite  element 
cracked  in  isotropic  tension  state  of  stress,  or  crushed  in  isotropic  compres- 
sion or  tension-compression  state  of  stress,  the  stiffness  contribution  of  that 
integration  point  to  the  element  stiffness  was  considered  as  zero.  If  all 
the  integration  points  of  all  the  elements  loaded  in  compression  were  crushed 
in  isotropic  compression  or  in  tension-compression,  the  failure  strength 
of  the  section  was  reached.  Thus,  depending  on  the  loading  on  the  section, 
as  a consequence,  the  depth  of  the  concrete  section  loaded  in  compression, 
the  failure  strength  of  the  section  would  vary.  Thus,  a series  of  values  of 
failure  strengths  for  the  section  could  be  obtained  as  a function  of  the 
loading  on  the  section,  resulting  in  an  envelope,  now  defined  as  the  failure 
strength  envelope  or  simply  the  failure  envelope. 
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The  main  purpose  of  the  analysis  then  was  to  obtain  a failure  envelope 


! 


for  each  type  of  section  analyzed.  In  order  to  achieve  this,  the  section 
was  imposed  with  monotonical ly  increasing  deformations  in  a fixed  ratio 
of  axial  to  rotational,  until  the  failure  strength  of  the  section  was 
attained. 

Since  the  ratio  of  imposed  deformations,  axial  to  rotational  was  fixed 
while  monotonical ly  increasing  the  value  of  these  deformations,  the  neutral 
axis  of  imposed  deformations  was  fixed.  Each  section  was  subjected  to  seven 
deformation  ratios  or  patterns:  pure  compression,  pure  extension,  pure 
rotation,  extension/rotation  of  5.0,  extension/rotation  of  15.0,  compression/ 
rotation  of  5.0  and  compression/rotation  of  15.0.  For  each  section,  seven 
failure  strengths  were  obtained  corresponding  to  the  seven  deformation  ratios 
or  patterns.  These  failure  strengths  were  plotted  to  give  the  failure 
envelope  for  the  section.  The  failure  envelope  was  plotted  in  terms  of 
deformations  which  were  reduced  to  dimensionless  quantities  of  strains  and 
curvatures.  Note  that  the  section  was  imposed  with  deformations  rather  than 
forces.  This  approach  was  adopted  since  the  final  aim  was  to  obtain  the 
failure  envelope  in  terms  of  deformations  rather  than  forces.  A typical 
failure  envelope  is  indicated  in  Figure  15.  As  the  deformations  on  the  section 
were  monotonically  increased,  at  each  load  step  increment,  the  axial  force 
and  bending  moment  on  the  section  were  also  calculated.  These  response 
quantities  were  plotted  as  moment-rotation  (M-e)  and  axial  force  deformation 
(P-6)  curves  for  each  section.  Note  that  the  M-0  and  P-6  curves  were  referenced 
to  the  geometric  centroid  of  the  section.  The  failure  envelopes  also  indicated 
strains  and  curvatures  at  failure  at  the  centroid  of  the  section.  Note  also 
that  only  confined  sections  were  studied. 
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2.  RECTANGULAR  SECTION  RESULTS 


The  total  number  of  cases  analyzed  for  the  confined  rectangular  section 
(CRS)  is  indicated  in  Table  1.  The  strength  of  concrete,  the  percentage  of 
longitudinal  reinforcement,  and  the  area  of  confining  steel  were  selected  as 
the  parameters  that  would  be  varied.  The  first  two  sections  analyzed  were 

CRS  #3  and  CRS  #7  (Table  1).  Each  of  the  above  two  sections  contained  2 percent 

longitudinal  reinforcement  both  in  tension  and  compression  zones  (p  = p1  = 2%). 
The  area  of  transverse  steel  per  leg  was  3.50  in2  over  a length  of  19.685  inches, 
calculated  on  the  basis  of  ACI  1971  code  requirements  for  seismic  resistant 
design.  The  only  difference  between  the  two  was  the  strength  of  concrete; 
the  ultimate  uniaxial  compressive  stress  of  concrete  for  CRS  #3  was  5 ksi , 

while  that  of  CRS  #7  was  6 ksi.  CRS  #4  and  #8  were  sections  with  p = p ' = 1/4%, 

while  the  ultimate  uniaxial  compressive  stresses  for  concrete  were  5 ksi  and 
6 ksi,  respectively.  CRS  #11  had  p = p ' =2%  and  an  ultimate  uniaxial  compres- 
sive stress  of  concrete  of  9.0  ksi.  CRS  #11.1  was  similar  to  CRS  #11  with 
the  exception  of  the  uniaxial  compression  stress-strain  relationship  for  the 
concrete,  which  was  different.  Section  CRS  #6  was  analyzed  with  1 percent 
longitudinal  reinforcement  at  each  face  and  an  ultimate  uniaxial  concrete 
compressive  stress  of  5.0  ksi.  While  keeping  the  section  size,  the  percentage 
of  longitudinal  reinforcement  and  the  properties  of  concrete  constant  and 
equal  to  that  of  CRS  #4,  the  transverse  reinforcement  was  varied  to  study  the 
effect  of  confinement  on  the  failure  envelope  (CRS  #4.2).  The  effect  of 
refinement  of  mesh  in  the  response  characteristics  and  the  failure  envelope 
was  studied  in  the  analysis  of  Section  CRS  #4.1. 
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TABLE  1 CONFINED  RECTANGULAR  SECTION  (CRS) 
PARAMETRIC  STUDY  LAYOUT 


a.  Analytical  Model inq  For  Concrete  And  Steel 

The  dimensions  of  the  section  were  kept  constant  throughout  the 
analysis.  No  cover  was  considered  in  the  analysis,  since  the  cover 
would  spall  at  a very  early  stage  of  loading,  thereby  not  contributing 
significantly  to  the  ultimate  capacity  of  the  section.  Inclusion  of 
the  cover  would  also  cause  numerical  stability  problems.  Concrete  was 
modeled  using  three-dimensional  node  finite  elements.  Most  of  the 
sections  were  analyzed  using  a four  layer  system.  Some  sections  were 
analyzed  using  eight  layers  of  brick  elements. 

Steel  was  modeled  as  three-dimensional  truss  elements.  These 
elements  were  connected  at  the  nodes  of  the  concrete  brick  elements. 

The  confined  section  under  consideration  is  shown  in  Figure  16.1  while 

an  analytical  model  is  shown  in  Figure  16.2. 

Boundary  conditions  imposed  on  the  section  were  as  follows: 
symmetry  boundary  conditions  were  imposed  on  the  face  1-11-15-5 
(Figure  16.2)  in  the  X-direction  while  nodes  10  and  20  were  restrained 
in  the  Y-direction.  To  prevent  rigid  body  modes  of  deformation  in  the 
Z-direction,  face  1-5-10-6  was  fixed  in  this  direction.  Nodes  3-8-18-13 
defined  the  horizontal  plane  passing  through  the  centroidal  axis  of  the 
section  (Figure  16.2).  Displacements  (extension  or  compression)  were 
imposed  on  the  plane  defined  by  nodes  6,  10,  16,  20,  while  this  same 
plane  (nodes  6-10-16-20)  was  subjected  to  rotations  about  the  horizontal 
plane  through  the  centroidal  axis  of  the  section  (nodes  3-8-18-13). 
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The  moment  and  axial  force  measured  as  response  quantities  were 
measured  at  the  geometric  centroidal  axis.  Thus,  the  M-e  curves  and  the 
P-6  curves  were  plotted  at  the  centroidal  axis  of  that  section. 

The  history  of  nonlinearities  that  occur  throughout  the  loading 
history  are  presented  as  events  in  tabular  form  on  the  plots  of  M-e 
and  P-6  curves.  Each  event  is  identified  by  a number  (for  example, 
refer  to  Figure  26.3).  The  events  are  described  at  the  gaussian 
integration  points  defined  in  these  figures  and  tables  as  "nodes". 

Thus,  if  a nonlinear  event  occurs  during  loading  (e.g.,  "cracking  of 
3 nodes  of  element  #4",  identified  as  point  1 in  Figure  26.3)  in  effect 
cracking  of  3 gaussian  integration  points  occurs.  Note  also  that  the 
current  status  of  each  element  is  defined.  Thus,  in  Figure  26.3,  at 
point  5 identified  in  this  figure,  7 gauss  points  of  element  #2  undergo 
yield  while  at  point  8 the  remaining  gauss  point  yields.  Thus  at  point  8, 
the  current  status  of  element  #2  is  that  all  8 gauss  points  of  element  #2 
have  yielded. 

(1)  Response  Of  CRS  #3 

Confined  rectangular  section  CRS  #3  had  concrete  with  an 
ultimate  uniaxial  compressive  stress  of  5.0  ksi . The  uniaxial 
stress  strain  curve  used  for  concrete  in  the  analysis  is  indicated 
in  Figure  17.  The  idealized  curve  is  shown  by  the  dashed  line. 

Note  that  concrete  behavior  is  assumed  elastic  and  linear  up  to 
the  point  where  the  stress  is  0.48  times  the  ultimate  compressive 
stress.  After  this  point  the  nonlinear  behavior  is  idealized 
by  three  linear  portions  of  nonelastic  behavior  as  shown  in 
Figure  17.  Ultimate  strain  was  assumed  to  be  0.003  in/in.  The 
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STRAIN  e (in/in) 

Idealized  model  used  in  analytical  response  prediction 

Typical  stress-strain  curve.  a 


a.  Obtained  from  Ref.  8,  pp.  12.  Ultimate  strain  was  assumed  to  be 
0.003  in/in. 


Note:  This  stress-strain  relation  used  for  all  cases  except 
CRS  #7,  CRS  #8,  CRS  #11,  CRS  #11-1 


Figure  17.  Stress-Strain  Curve  for  Concrete 
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Figure  18  Stress-Strain  Curve  for  Steel.  Grade  60  ASTM 


stress-strain  curve  used  for  steel  is  shown  in  Figure  18.  It 
was  idealized  as  a bilinear  curve  with  a strain  hardening  ratio 
of  0.02;  i.e.,  Es^/Es  = 0.02  where  Esh  is  the  slope  of  the  strain 
hardening  curve  and  Es  is  the  Initial  modulus  of  elasticity  of 
steel.  The  failure  envelope  for  this  section  using  the  above 
material  properties  is  shown  in  Figure  23.  while  the  response  of 
this  section  to  various  deformation  ratios  is  indicated  in 
Figures  24.1  through  24.7  in  terms  of  M-e  and  P-6  curves.  Inter- 
pretation of  each  of  these  response  curves  is  discussed  for  this 
section  in  detail.  For  concrete  sections  presented  later  in 
this  section  of  the  report  only  the  main  differences  are  discussed. 

Figure  24.1  shows  the  plot  of  P vs.  6 for  a pure  extension  case. 
The  response  is  linear  up  to  the  point  prior  to  1 indicated  in 
Figure  24.1.  At  point  1 cracking  of  the  whole  section  occurs  and 
the  total  resistance  is  then  due  to  the  longitudinal  steel.  At 
point  2 in  Figure  24.1,  the  yield  of  longitudinal  steel  occurs. 
Failure  of  this  section  under  pure  extension  is  then  controlled  by 
serviceability  rather  than  strength.  A ductility  ratio,  defined 
as  a ratio  of  maximum  strain  in  longitudinal  steel  to  the  initial 
yield  strain  of  longitudinal  steel,  of  10  to  15  would  be  an  appro- 
priate value  where  failure  would  be  defined.  However,  on  the  plot 
of  failure  envelope  for  this  section  (Figure  23),  this  point  is 
not  plotted. 

In  Figure  24.2,  the  response  P versus  the  deformation  6 for  a 
case  of  imposed  pure  compressive  deformation  is  indicated.  Here  the 
response  is  linear  up  to  the  point  1 where  the  yield  of  concrete 


section  occurs  in  compression.  Yield  of  the  section  is  defined 
when  the  effective  stress  ag  (equation  1)  is  greater  than  t0,  the 
initial  effective  yield  stress  (equation  7).  At  point  2 in 
Figure  24.2,  the  deformations  have  been  increased  to  a point  where 
the  yield  of  longitudinal  bars  occurs  in  compression.  In  the  next 
load  step,  after  the  yield  of  re-bars  has  occurred,  crushing  of 
concrete  occurs  in  compression.  Failure  of  t^e  section  was  defined 
at  this  point,  which  is  indicated  on  the  failure  envelope 
(Figure  23). 

The  response  of  the  section  to  pure  rotation  imposed  deformations 
is  shown  in  Figure  24.3.  The  section  remains  elastic  up  to  point  1. 

In  the  load  step  defined  by  point  2,  cracking  of  the  bottom  element 
occurs.  As  deformations  are  increased  further,  the  bottom  half  of 
the  section  cracks  completely  in  tension  [element  nos.  3,  and  4, 
(Figure  16.2)].  Note  that  the  cracking  of  elements  is  followed  by 
a loss  of  stiffness  as  indicated  by  the  M-e  curve.  As  deformations 
are  increased  further,  certain  integration  points  of  the  top 
element  #1  crush  in  tension-compression.  Following  the  yield  of 
re-bars  in  tension  and  compression  (point  8 in  Figure  24.3), 
successive  crushing  of  gauss  points  of  elements  1 and  2 occurs  in 
isotropic  compression  and  in  tension-compression  until  at  point  13, 
all  the  elements,  loaded  in  compression,  crush  completely.  After 
this  point  the  only  resisting  force  is  due  to  the  strain  hardening 
effect  of  steel.  Failure  of  the  section  for  this  imposed  deformation 
is  thus  defined  when  all  the  elements  of  concrete  loaded  in  compres- 
sion crush  completely  in  an  isotropic  compression  or  a tension- 
compression  state  of  stress.  The  deformation  (curvature)  associated 
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with  the  section  at  this  point  is  shown  on  the  failure  envelope  in 
Figure  23. 

For  the  imposed  deformation  ratio  of  extension/rotation  of  5.0, 
the  deformations  imposed  were  such  that  the  neutral  axis  lies  between 
the  element  boundaries  of  element  #2  (Figure  19).  Thus,  element  #1 
and  top  portion  of  element  #2  are  in  the  compression  zone  while  the 
bottom  portion  of  element  #2,  as  well  as  all  of  elements  #3  and 
#4  are  in  tensile  zone.  The  position  of  this  neutral  axis  is  now 
fixed  while  the  deformations  are  monotonically  increased.  It  can  be 
seen  from  Figure  24.4  that  the  axial  force  on  the  element  is 
initially  tensile  in  nature.  Increase  of  the  imposed  deformations 
results  in  the  cracking  of  elements  #3  and  #4  and  the  bottom  four 
integration  points  of  element  #2  (point  3 in  Figure  24.4).  However, 
in  the  next  load  step,  the  top  four  integration  points  yield  in 
tension-compression,  an  event  which  causes  an  out-of-balance  force 
in  the  solution  process.  The  process  of  iterating  the  solution  to 
convergence  causes  the  stress  state  to  change  to  a pure  tension  state, 
resulting  in  the  cracking  of  the  top  four  nodes  of  the  element  #2. 

Once  the  cracking  of  the  tensile  zone  has  occurred,  most  of  the 
forces  resulting  from  the  imposed  deformations  are  resisted  by 
the  longitudinal  tensile  steel,  which  is  still  elastic,  while  a 
small  portion  is  resisted  by  compressive  steel  and  the  concrete  in 
the  compressive  zone.  Hence,  the  axial  force  in  the  section  is  still 
tensile  and  increases  with  the  increase  in  deformation.  As  soon  as 
the  tensile  steel  yields  in  tension,  the  forces  resulting  from  a 
further  increase  in  deformation  are  resisted  mainly  by  the  steel  in 
compression  and  the  concrete  compressive  block.  The  increase  in 
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force  in  the  tensile  steel  is  very  small.  Such  a force  distri- 
bution causes  a drop  in  the  axial  force  which  is  tensile  and 
eventually  leads  to  a compressive  axial  force  P on  the  section. 

At  point  10  in  Figure  24.4,  the  re-bars  in  compression  yield. 
Further  increases  in  deformation  cause  crushing  of  the  integration 
points  of  the  concrete  in  the  compressive  zone  until  at  point  13, 
crushing  of  the  whole  concrete  compressive  block  occurs.  At  this 
stage  any  increase  in  deformation  is  resisted  by  the  steel  in 
tension  and  compression,  both  of  which  have  already  undergone 
yield.  Failure  of  the  section  for  this  deformation  ratio  is 
defined  by  point  13,  where  all  the  concrete  loaded  in  compression 
has  crushed,  which  is  shown  on  the  plot  of  the  failure  envelope 
in  terms  of  strains  and  curvatures  in  Figure  23. 

Response  of  the  section  when  subjected  to  a monotonically 
increasing  deformation  ratio  of  extension/rotation  of  15  is  indi- 
cated in  Figure  24.5.  In  this  case,  the  neutral  axis  of  deforma- 
tions lies  outside  the  section,  thereby  forcing  the  whole  section 
to  be  in  tension  (Figure  20).  Thus,  the  case  is  a trivial  case 
where  the  increase  in  deformation  eventually  cracks  the  whole 
section  in  isotropic  tension  (point  3 in  Figure  24.5).  After 
the  section  has  cracked,  any  additional  deformation  is  resisted 
by  the  longitudinal  steel  which  is  elastic  up  to  point  4,  where 
the  bottom  re-bars  yield  in  tension.  Hence,  the  response  is 
controlled  by  steel  rather  than  concrete,  and  failure  in  this 
case  is  defined  similar  to  failure  in  the  pure  extension  case. 

The  failure  point  has  not  been  plotted  for  the  deformation  ratio 
on  the  failure  envelope  in  Figure  23. 
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Deformation  ratio  of  compression/rotation  of  5.0  is  opposite 
to  extension  rotation  of  5.0  as  far  as  the  location  of  the  neutral 
axis  of  deformation  is  concerned.  Element  #4  and  the  bottom  por- 
tion of  element  #3  are  in  tension,  while  elements  #1,  #2  and  the 
top  portion  of  element  #3  are  in  compression  (Figure  21). 

The  response  of  the  section  to  this  deformation  pattern  is 
indicated  in  Figure  24.6.  As  a contrast  to  the  response  of  the 
section  to  a deformation  pattern  of  extension/rotation  of  5,  it 
can  be  observed  that  the  axial  force  on  the  member  is  always 
compressive  in  nature.  This  behavior  is  expected  as  the  compressive 
zone  is  the  main  contributor  to  the  resisting  mechanism.  Once  the 
bottom  element  has  cracked  in  tension,  causing  some  loss  of  initial 
stiffness,  future  behavior  to  increase  in  deformation  is  controlled 
by  the  steel  and  the  concrete  in  the  compressive  zone.  With  suc- 
cessive increase  in  deformation  the  concrete  in  the  compressive 
zone  crushes.  Failure  of  the  section  to  this  deformation  pattern 
is  defined  when  the  concrete  section  in  the  compressive  zone  is 
completely  crushed.  Such  a deformation  condition  is  plotted  on 
the  failure  envelope  of  Figure  23. 

The  last  deformation  pattern  imposed  on  the  section  was  a 
combination  of  compressive  displacement  to  rotation  of  15.  This 
deformation  pattern  was  such  that  the  neutral  axis  of  imposed 
deformations  was  outside  the  section,  thereby  imposing  compressive 
deformations  on  the  section  (Figure  22)-  The  response  of  the  sec- 
tion is  indicated  in  Figure  24.7.  As  can  be  seen  from  this  figure, 
as  the  deformations  are  monotonically  increased,  certain  elements 
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Figure  23  CRS  #3  Failure  Envelope 


Fig.  24.1  CRS  #3  - Pure  Extension 
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yield  in  isotropic  compression,  while  some  yield  in  tension- 
compression.  Further  successive  increases  in  deformation  cause 
the  top  re-bars  to  yield  in  compression  and  then  the  bottom  re-bars 
to  yield  in  compression.  Failure  for  this  section  is  defined  when 
the  concrete  section,  loaded  in  compression,  totally  crushes  either 
in  isotropic  compression  or  tension-compression  state  of  stress. 

The  deformations  associated  with  this  failure  point  are  indicated 
in  Figure  23. 


(2)  Response  Of  CRS  #4 

CRS  #4  has  the  same  stress-strain  relationship  for  concrete  and 

steel  properties  as  CRS  #3  (Figures  17  and  18).  Four  layers  of 

bricks  were  used  to  model  the  concrete.  Longitudinal  reinforcement 

was  1/4%  at  each  face.  Since  the  transverse  reinforcement  required 

by  AC  I 1971  code  for  seismic  design  requirements  is  independent  of 

the  percentage  of  longitudinal  reinforcement,  but  is  dependent  on 

the  size  of  the  section,  spacing  between  stirrups,  ultimate  uniaxial 

compressive  stress  of  concrete  f*  and  the  yield  stress  of  steel  f , 

c y 

the  area  of  transverse  reinforcement  to  be  provided  for  this  section 
would  have  been  the  same  as  for  section  CRS  #3.  This  would  have 
created  a situation  where  the  area  of  transverse  reinforcement  would 
be  greater  than  the  area  of  longitudinal  reinforcement.  Such  a situ- 
ation would  be  an  unrealistic  one  and  the  problem  would  be  only 
academic  in  nature.  Hence,  it  was  decided  instead  to  use  an 

i 

arbitrary  area  of  transverse  reinforcement  which  would  seem  reason- 
able compared  to  the  area  of  longitudinal  reinforcement.  Since 
the  area  of  confining  steel  is  so  small,  as  is  also  the  area  of 
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longitudinal  steel,  this  section  is  similar  to  an  unconfined 
section.  The  response  of  the  section  to  imposed  deformations 
would  be  controlled  by  the  behavior  of  concrete,  with  little 
effect  of  the  behavior  of  longitudinal  steel.  In  effect,  it  is 
very  close  to  a plain  concrete  model.  Since  the  confining  effect 
of  transverse  reinforcement  is  almost  negligible,  the  value  of 
Poisson's  ratio  was  set  to  10"7.  This  also  stabilized  the 
numerical  solution  process.  The  failure  envelope  for  this 
section  is  indicated  in  Figure  25,  while  the  response  of  the 
section  to  various  deformation  ratios  is  presented  in  Figures 
26.1  through  26.7. 

In  Figure  26.1,  the  pure  extension  imposed  displacements  case, 
it  can  be  seen  that  the  initial  stiffness  and  load  carrying  capacity 
is  mainly  due  to  the  concrete  section.  Once  the  section  completely 
cracks,  the  resistance  to  imposed  deformation  is  due  to  steel. 
Comparing  the  CRS  #4  response  to  that  of  Figure  24.1  of  section 
CRS  #3,  Figure  24.1  reveals  that  in  the  response  of  CRS  *3, 
longitudinal  steel  was  the  main  source  of  strength,  and  the 
resisting  force  was  mainly  due  to  the  forces  in  steel. 

In  Figure  26.2,  where  the  imposed  deformations  are  pure 
compressive  displacements,  it  can  be  observed  that  the  crushing 
of  concrete  in  isotropic  compression  occurs  before  the  yield  of 
longitudinal  bars,  unlike  the  failure  of  CRS  #3  when  subjected 
to  a similar  deformation  pattern.  Failure  of  the  section  in  this 
case  was  also  defined  by  point  2 (Figure  26.2)  where  the  concrete 
section  totally  crushes  in  compression-compression. 
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Figure  26.3,  the  pure  rotation  imposed  deformations  case, 
indicates  all  the  events  that  occur  throughout  the  deformation 
history.  Failure  of  the  section  was  defined  when  all  the  concrete 
loaded  in  compression  crushed  in  isotropic  compression  or  tension- 
compression  stress  states.  Note  that  the  longitudinal  steel  yields 
before  crushing  of  the  concrete  section  occurs. 

Failure  of  the  section  when  subjected  to  an  imposed  deformation 
pattern  of  extension/rotation  = 5.0  occurs  before  the  longitudinal 
steel  yields.  Failure  is  defined  by  the  deformation  at  which  the 
concrete  loaded  in  the  compression  zone  totally  crushes  (Figure  26.4). 

In  Figure  26.7,  which  shows  the  response  of  the  section  to 
an  imposed  deformation  pattern  of  compression/rotation  of  15,  failure 
of  the  concrete  section  occurs  after  the  yield  of  top  re-bars  in 
compression.  Note  that,  unlike  the  response  of  CRS  #3  to  a similar 
deformation  ratio  (Figure  24.7),  the  bottom  re-bars  of  Section  CRS 
#4  do  not  yield. 

The  failure  envelope  of  Section  CRS  #4  shown  in  Figure  25  is 
quite  different  from  the  one  indicated  in  Figure  19  for  Section 
CRS  #3.  Section  CRS  #4  behaves  essentially  like  an  unconfined  plain 
concrete  section.  Since  behavior  is  mainly  controlled  by  the 
behavior  of  concrete,  the  response  is  very  sensitive  to  the  mesh 
size.  A four-layer  mesh  may  be  a coarse  mesh  which  does  not 
accurately  account  for  the  nonlinear  behavior  through  the  depth  of 
the  section.  This  could  be  the  cause  of  the  difference  in  shape  of 
the  failure  envelopes  of  Sections  CRS  #3  and  CRS  #4.  Hence,  section 
CRS  #4  was  reanalyzed,  refining  the  mesh,  using  an  eight  layer  mesh 
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Figure  26.1  CRS  #4-  Pure  Extension 
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Instead  of  a four  layer  mesh.  Results  of  this  analysis  are 
presented  next. 

(3)  Response  of  CRS  #4.1 

CRS  4.1  had  similar  properties  of  steel  and  concrete  and 
the  same  amount  of  longitudinal  and  transverse  reinforcement  as 
CRS  #4.  An  eight  layer  mesh  was  used  to  model  the  concrete  as 
shown  In  Figure  27.  The  failure  envelope  (Figure  28)  is  a great 
improvement  over  that  of  Figure  25  for  CRS  *4.  Response  of 
the  section  to  various  deformation  patterns  Is  Indicated  in 
Figures  29.1  through  29.6. 

As  can  be  seen  from  Figure  29.1,  for  pure  compression  imposed 
displacements,  a more  distinct  and  refined  P-6  curve  is  obtained 
with  an  eight  layer  mesh,  with  a distinct  negative  slope  after  the 
maximum  load  point  is  reached. 

The  most  significant  change  in  response  is  in  the  case  where 
the  section  is  imposed  with  pure  rotation  deformations  as  indicated 
in  Figure  29.2  compared  to  Figure  26.3.  Failure  in  this  case 
occurs  at  a lower  magnitude  of  rotation  of  0.0017  radians. 

The  deformation  at  which  failure  of  the  section  occurs  when 
subjected  to  extension/rotation  of  5.0  is  close  to  the  one  predicted 
with  a four  layer  mesh  (Figure  29.3  compared  to  Figure  26.4).  The 
response  to  a deformation  pattern  of  extension/rotation  of  15  is 
also  similar  (Figure  29.4  compared  to  Figure  26.5). 
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Figure  27  Eight  Layer  Finite  Element  Model 
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However,  the  response  and  the  deformation  at  failure  for  an 
imposed  deformation  of  compress ion/ rotation  of  5 is  different. 

As  can  be  seen  from  Figure  29.5,  the  eight  layer  model  predicts 
more  strength  than  the  four  layer  model  (Figure  26.6).  The  eight 
layer  model  is  numerically  a more  stable  system.  The  nonlinearities 
and  consequent  changes  in  stiffness  are  not  as  drastic  as  in  the 
four  layer  system. 

Figure  29.6,  however,  indicates  that  the  failure  strength  of 
the  section  is  smaller  when  subjected  to  a deformation  pattern  of 
compression/rotation  of  15.0.  This  could  be  due  to  the  multi- 
axial  stress  behavior  where  a larger  portion  of  the  section  fails 
in  tension-compression  In  the  eight  layer  model  as  compared  to 
the  four  layer  model. 

It  should,  however,  be  emphasized  at  this  stage  that  whenever 
the  response  of  the  section  is  predominantly  controlled  by  the 
behavior  of  concrete,  a finer  mesh  would  give  more  accurate  results 
and  a more  stable  solution  process. 

(4)  Response  Of  CRS  #6 

CRS  #6  was  analyzed  with  a longitudinal  reinforcement  of  1% 
on  each  face  (p  = p'  = 1%)  and  with  concrete  having  an  ultimate 
compressive  stress  of  5.0  ksi.  The  amount  of  transverse  reinforcement 
was  assumed  to  be  40%  of  the  longitudinal  reinforcement.  A section 
so  reinforced  would  be  expected  to  have  a high  ductility. 

The  failure  envelope  for  CRS  #6  is  shown  in  Figure  30  iwhile 
the  response  in  terms  of  M-0  and  P-6  curves  is  shown  in  Figures  31.1 
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through  31.6.  This  failure  envelope  shows  that  the  failure  strength 
of  the  section  CRS  #6  is  higher  when  subjected  to  pure  compressive 
displacements  as  compared  to  CRS  #3  with  a reinforcement  ratio  of 
2%  (Figure  23  ),  when  subjected  to  similar  compressive  displacements. 
Figure  30  (compared  to  Figure  23)  also  shows  that  the  failure 
strength  of  the  section  is  higher  when  subjected  to  deformation 
ratios  of  compression/rotation  of  5 and  compression/rotation  of  15, 
while  the  failure  strength  is  decreased  when  subjected  to  a deforma- 
ratio  of  extension/rotation  of  5.0.  The  deformation  capacity  of 
the  section  remains  the  same  when  subjected  to  a pure  rotation. 

Note  that  the  value  of  Poisson's  ratio  for  CRS  #6  was  10"6  as  compared 
to  a value  of  0.2  for  CRS  #3. 

The  general  pattern  of  the  M-9  and  P-6  curves  in  Figures  31.1 
through  Figure  31.6  is  the  same  as  that  of  CRS  #3  when  subjected  to 
corresponding  deformation  ratios.  Thus,  in  the  first  stage  of 
study  the  uniaxial  compressive  stress-strain  relationship  was  kept 
the  same  while  the  effect  of  different  parameters  were  studied; 
e.g.,  change  in  longitudinal  reinforcement  (CRS  #4,  CRS  #6),  change 
in  mesh  size  (CRS  #4,  CRS  #4.1),  and  change  in  transverse  reinforcement 
(CRS  #4,  CRS  #4.2,  the  response  of  which  is  discussed  later  In  this 
section). 


Hence,  higher  strength  concrete  was  used,  and  only  the  amount 
of  longitudinal  steel  was  varied  to  study  the  effect  on  the  shape 
of  the  failure  envelope.  Two  different  strengths  of  concrete  were 
used;  one  with  an  ultimate  uniaxial  compressive  stress  of  6.0  ksl 

and  the  other  with  an  ultimate  uniaxial  compressive  stress  of 
9.0  ksl.  Results  of  these  analyses  are  presented  next. 
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Figure  31. 3 CRS  #6  - Pure  RoUttoo 


EVENT 
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Figure  31.4  CRS  #6  - Extenslon/Rotetlon  * 5.0 


Flgvrt  31.7  CHS  #4  - Compression/ Rotation  • lf.0 


(5)  Response  Of  CRS  #7 


CRS  #7  was  similar  to  CRS  #3  in  dimensions,  percentage  of 
longitudinal  reinforcement,  area  of  transverse  reinforcement, 
and  the  properties  of  steel  (stress-strain  relation).  The  ily 
difference  between  the  two  sections  was  the  uniaxial  stress-strain 
curve  for  concrete.  Figure  32  Shows  the  stress-strain  curve 
used  in  the  analysis  of  CRS  #7. 

The  failure  envelope  plotted  from  the  results  of  the  analyses 
of  the  section  subjected  to  various  deformation  ratios  is  shown 
in  Figure  33.  This  envelope  is  very  similar  in  shape  to  that  of 
CRS  *2  (Figure  23). 

Figures  34.1  through  34.6  indicate  that  the  M-0  and  P-6  curves 
for  CRS  #7  are  similar  to  those  of  CRS  #3,  and  the  modes  of  failure 
are  also  similar. 


(6)  Response  of  CRS  #8 

Section  CRS  #8  was  similar  to  CRS  #4,  in  dimensions,  amount 
of  longitudinal  reinforcement,  and  amount  of  transverse  reinforce 
ment.  The  stress-strain  curve  used  to  represent  the  properties 
of  concrete  is  shown  in  Figure  32. 


Results  of  the  analyses  of  this  section  subjected  to  various 
deformation  ratios  are  plotted  in  the  form  of  a failure  envelope 
in  Figure  35.  Comparing  this  envelope  to  that  of  CRS  #4  (Figure  25) 
shows  a distinct  difference  between  the  two  curves,  especially  for 
the  case  when  the  section  is  subjected  to  pure  rotation  deformation. 


STRESS  (ksi) 


Actual  Curve0 

Idealized  Curve  for 
Analysis 


1 r—  I 

0.002 

STRAIN  (in/in) 


a)  Curve  obtained  from  ref.  8,  p.  12.  Ultimate  strain  was 
assumed  to  be  0.003  in/in. 

Ref.  8:  Park,  R. , and  Paulay,  T. , Reinforced  Concrete  Structures 
John  Wiley  & Sons,  Inc. , 1975 


Figure  32  Stress-Strain  Curve  for  Concrete.  Ultimate  Stress  = 60  ks 
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Compression 


Fig.  33  CRS  *1  Failure  Envelope 


Compression  Displacement  « (10"*  Ins) 


Figure  34.2  CHS  #7  - Pure  compression 


The  failure  strength  estimated  by  the  analysis  of  CRS  #4  is 
significantly  higher  than  that  of  CRS  #8  (Figure  35).  A 
comparison  with  that  of  CRS  #4.1  (Figure  38)  shows  a better 
resemblance  between  the  two  failure  envelopes.  The  shape  and 
magnitude  of  the  two  curves  follow  a close  pattern. 

Figures  36.1  through  36.7  indicate  the  M-0  and  P-6  curves 
for  this  section.  These  curves  also  follow  a close  resemblance 
to  those  of  CRS  #4.1  (Figures  29.1  through  29.6).  Modes  of 
failure  are  also  similar  for  various  deformation  patterns  in 
both  cases. 

(7)  Response  Of  CRS  #11 

It  was  decided  to  study  the  effect  of  the  strength  of 
concrete  on  the  failure  envelope.  As  an  upper  bound,  concrete 
with  an  ultimate  uniaxial  compressive  stress  of  9.0  ksi  was  used 
with  a longitudinal  steel  ratio  of  2 % at  each  face.  The  transverse 
reinforcement  area  used  was  3.50  in2.  Since  no  experimental  curves 
were  available  for  such  a high  strength  concrete,  an  analytical 
curve  was  used.  Figure  37  shows  the  stress-strain  curve  as 
suggested  by  Saenz  In  the  paper  by  Popovics  (Ref.  7). 

Figure  38  shows  the  failure  curve  of  such  an  analysis.  The 
figure  indicates  a capacity  of  the  section  when  subjected  to  pure 
rotation  to  be  less  than  that  of  CRS  #3  with  a lower  strength  concrete 
(Figure  23).  The  capacities  of  the  section  when  subjected  to 


7.  Popovics,  S.,  "A  Review  of  Stress-Strain  Relationships  for  Concrete", 
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Figure  34.3  CRS  47  - Extension/Rotation 
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COMPRESSIVE  DISPLACEMENT  i (10‘*  Ins) 


Figure  34.6  CRS  #7  - Compression/Rotation 


Compression 


Figure  35  CRS  #8  - Failure  Envelope 


EVENT 


EXTENSION  6 


COMPRESSION 


Figure  36.2  CRS  #8  - Pure  Compression 


Figure  36.3  CDS  IS  - Pure  PoUtlon 


Figure  36.4  CRS  #8  - Extenilon/RoUtlon  • 5.0 
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36.5  CRS  #8  - Extension/Rotation  * 15.0 


Figure  36.6  CHS  #8  - Compress  Ion/ Rotation  ■ 5.0 


#8  - Coapress tofl/RoUtton  » )5.0 


deformation  ratios  of  extension/rotation  of  5,  compression/ 
rotation  of  5,  compression/rotation  of  15,  and  pure  compression 
are  higher  than  that  of  CRS  #3  (Figure  19).  Possibly  a more 
refined  eight  layer  mesh  would  result  in  a better  estimation  of 
the  capacity  when  subjected  to  pure  rotation. 

Figures  39.1  through  39.7  shows  the  response  curves  (M-e 
and  P-6  curves)  for  various  deformation  ratios.  Again,  the  general 
form  of  these  curves  is  similar  to  CRS  #3. 

(8)  Response  Of  CRS  #11.1 

CRS  #11.1  was  a section  with  uniaxial  compressive  strength  of 
9.3  ksi,  2%  longitudinal  reinforcement  at  each  face  and  a transverse 
reinforcement  area  ASh  of  3.50  ins2.  This  section  was  analyzed  using 
a refined  eight  layer  finite  element  mesh.  The  uniaxial  stress- 
strain  curve  used  for  concrete  is  shown  in  Figure  40.  It  should 
be  noted  that  the  ultimate  uniaxial  concrete  compressive  strain  used 
is  0.0045  in/in,  as  compared  to  the  ultimate  uniaxial  concrete 
compressive  strain  of  0.003  in/in  used  to  analyze  all  the  other 
rectangular  sections. 

Failure  envelope  obtained  from  such  an  analysis  is  shown  in 
Figure  41.  It  can  be  seen  that  the  general  shape  of  this  failure 
envelope  in  the  compression  zone  is  similar  to  that  of  the  other 
rectangular  section.  Note  the  increase  in  deformation  strength  of 
CRS  #11.1  in  the  compression  zone  (when  subjected  to  deformation 
ratios  of  pure  compression,  compression/rotation  of  5.0,  pure 
rotation)  as  compared  to  similar  deformation  strength  of  sections 
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STRESS  (ksi ) 


0.0006  0.001  0.002  0.003 


STRAIN  (In/in) 

Ref.  7:  Popovics,  S.,  “A  Review  of  Stress-Strain  Relationships 
for  Concrete,"  Journal,  American  Concrete  Institute, 
March,  1970. 


Figure  37  Stress-Strain  Curve  for  Concrete.  Ultimate  Stress  * 9.0  ksi 
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COMPRESSION  DISPLACEMENT  6 (lO-4  Ins) 


AXIAL  LOAD  P,  (COMPRESSION)(k1ps)  MOMENT  M (klp-ln.) 
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Elems  #3  & 4 
cracked  completely. 

© 

Elem  #2  yield  4 nodes. 

© 

Elem  #2  yields  3 nodes. 
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Elem  #1  crushes  2 nodes. 

Elem  #1  crushes  4 nodes. 

NOTE 

NO. 

EVENT 

© 

Elem  #2  crushes  5 nodes. 

Elem  #2  cracks  2 nodes. 
Solution  does  not  converne. 

© 

Elem  *1  crushes  5 nodes. 
Steel  yields. 

© 

Elem  *1  crushes  5 nodes. 

© 

Elem  #1  crushes  completely. 
Solution  does  not  converge. 
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Figure  39.3  CRS  #11  - Pure  Rotation 
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Figure  39.5  CRS  #11  - Extension/Rotation  * 15.0 


Figure  39.6  CHS  #11  - Coopreiil on/Rotet Ion  * 5.0 
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with  lower  strength  concrete  (CRS  #3,  Figure  23;  CRS  #7,  Figure  33). 

This  increased  deformation  strength  is  due  to  the  higher  uniaxial 
ultimate  compressive  strain  of  0.0045  in/in  used  to  analyze  CRS 
#11.1.  Figure  42.1  through  42.6  indicate  the  response  of  CRS 
#11.1  when  subjected  to  various  deformation  ratios. 

(9)  Response  Of  CRS  #4.2 

CRS  #4  (confined  rectangular  section,  f'c  = 5.0  ksi,  p = p'  = 1/4*) 
was  reanalyzed  using  a higher  area  of  transverse  reinforcement.  The 
area  of  transverse  reinforcement  used  in  analysis  was  3.50  in2. 

Failure  envelope  of  such  a section  is  shown  in  Figure  43  while 
Figure  44.1  through  Figure  44.7  indicate  the  P-6  and  M-e  curves  for 
various  deformation  ratios.  The  concrete  section  was  modeled  as 
3-D  brick  elements  (four  layer  mesh.  Figure  16)  to  analyze  the 
response  for  all  deformation  ratios  except  pure  rotation.  Response 
of  the  section  when  subjected  to  pure  rotation  was  analyzed  using 
an  eight  layer  mesh. 

Comparing  the  failure  envelopes  of  CRS  #4  (transverse  reinforce- 
ment = 0.60  in2)  to  that  of  CRS  #4.2  (transverse  reinforcement 
= 3.50  in2),  one  can  observe  a higher  ductility  of  CRS  #4.2,  a 
phenomenon  to  be  expected  with  increase  in  transverse  reinforcement. 

It  should  be  mentioned,  that  the  problem  using  a four  layer  mesh  is 
numberically  very  sensitive.  Analysis  with  an  eight  layer  mesh 
might  have  resulted  in  a smoother  curve. 
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Compression 


CURVATURE  (TO-5  radians/in) 
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COMPRESSIVE  DISPLACEMENTS  (Ins) 
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Figure  43.  CRS  #4.2  Failure  Envelope 


Figure  44.1  CRS  #4.2  Pure  Extension 


COMPRESSIVE  DISPLACEMENTS  (ID-*  Ins)  NOTE 


Figure  44.2  CRS  #4.2  - Pure  Compression 
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Figure  44.5  CRS  #4.2  - Extension/Rotation 
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Figure  44.7  CRS  14.2  - Compress 1on/Rot4t1on  • 15.0 


3.  SUMMARY  OF  RECTANGULAR  SECTION  RESULTS 


The  response  of  various  sections  subjected  to  different  types  of 
loading  patterns  has  been  presented  in  detail  in  the  preceeding  para- 
graphs. Some  generalized  observations  are  now  summarized. 

(1)  The  shape  of  the  failure  envelopes  was  similar  for  all  the 
rectangular  sections  analyzed. 

(2)  The  numerical  solution  was  quite  sensitive  for  pure  rotation 
cases  and  lightly  reinforced  sections  with  transverse  steel.  Step  sizes  had 

to  be  reduced  and  finer  meshes  employed.  In  spite  of  the  extra  cost  we 
would  recommend  for  future  studies  that  finer  meshes  be  employed  for  all 
cases.  More  consistent  and  smoother  curves  would  result. 

(3)  When  the  strength  of  concrete  in  uniaxial  compression  was  increased 
from  5.0  ksi  to  6.0  ksl,  while  the  longitudinal  and  transverse  reinforcement 
was  kept  constant,  the  deformation  strength  of  the  section  reduced  with  an 
increase  of  concrete  strength.  (CRS  #3,  CRS  #7) 

(4)  The  deformation  strength  of  the  section  reduced  when  the  percentage 
of  longitudinal  steel  was  reduced  from  2 percent  to  1/4  percent.  Both  sections 
had  the  same  strength  concrete.  (CRS  #3,  CRS  #4.1) 

(5)  An  increase  in  the  amount  of  transverse  reinforcement  increased  the 
deformation  capacity  and  the  ductility  of  the  section.  (CRS  #4.1,  CRS  #4.2) 
Note  that  the  percentage  of  longitudinal  reinforcement,  and  the  strength  of 
concrete  in  uniaxial  compression  was  the  same  for  both  these  sections. 

4.  CYLINDRICAL  SECTION  RESULTS 

Development  of  failure  envelopes  for  cylindrical  sections  similar  to 
those  of  the  rectangular  sections  formed  the  second  major  phase  of  the  study 


of  failure  criteria  for  reinforced  concrete  structures.  Only  confined 
sections  were  analyzed.  The  main  parameter  of  interest  was  the  r/t  ratio 
of  the  cylindrical  section  and  its  effect  on  the  failure  envelope.  Two 
r/t  ratios  were  analyzed:  one  r/t  of  4 and  the  other  r/t  of  8.  (Note 
r = inside  radius  of  the  cylinder  and  t = thickness  of  the  cylinder). 

The  percentage  of  total  longitudinal  reinforcement  at  any  section  was 
2 percent.  Ultimate  uniaxial  compressive  stress  of  concrete  used  was 
5.0  ksi  (Figure  17)  while  the  properties  of  steel  were  idealized  as  shown 
in  Figure  18.  Both  cylindrical  sections  were  analyzed  for  seven  deformation 
ratios,  similar  to  the  rectangular  sections.  A third  cylindrical  section 
was  also  analyzed  where  the  section  with  r/t  ratio  of  4 was  subjected  to  a 
pressure  load  of  3-4-5.  The  cases  analyzed  are  shown  in  Table  2. 

a.  Analytical  Modeling  For  Concrete 

The  dimensions  of  the  cylindrical  section  analyzed  are  shown  in 

Figure  45.1.  The  thickness  of  the  section  was  constant  (1  meter  or  39.37 
inches)  while  the  radius  of  section  was  changed  depending  on  the  r/t 
ratio.  Thus  section  CCS- 1 had  a radius  of  157.48  inches  (4  meters)  while 
section  CCS-2  had  a radius  of  314.96  inches  (8  meters).  Since  the  section 
was  cylindrical,  and  the  loading  was  symmetric  about  its  diameter,  for 
the  purposes  of  analysis,  half  of  the  section  was  considered.  The  finite 
element  model  used  in  the  analysis  is  shown  in  Figure  45.2.  The  left  face 
of  the  section  was  fixed  while  deformations  were  imposed  on  the  right 
face  (face  24-28-32-36-40-37-33-29-25-21)  in  Figure  40. 
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(314.96") 


Table  2 

CONFINED  CYLINDRICAL  SECTION  (CCS) 


(1)  Response  Of  CCS-1 


Figure  46  shows  the  failure  envelope  of  confined  cylindrical 
section  CCS-1.  As  can  be  seen  from  the  figure,  the  failure  points 
for  deformation  ratios  of  extension/rotation  of  5 and  15  and  com- 
pression/rotation of  5 and  15  are  all  clustered  in  a single  area. 
This  is  because  the  neutral  axis  of  imposed  deformation  always  lies 
within  the  section  and  within  the  first  set  of  elements  close  to  the 
geometric  neutral  axis.  Again,  the  failure  strength  of  the  section 
is  less  when  subjected  to  pure  rotation  than  when  subjected  to 
deformation  ratios  of  displacement/rotation  of  5.0. 


Figures  47.1  through  47.6  indicate  the  response  in  terms  of 
M-0  and  P-6  curves  for  various  deformation  ratios.  A complete 
history  of  events  Is  also  indicated  throughout  the  loading  history 
for  each  deformation  ratio. 

(2)  Response  Of  CCS-2 

The  failure  envelope  for  CCS-2  with  r/t  ratio  of  8 is  indicated 
in  Figure  48.  Comparing  this  envelope  to  that  of  CCS-1  (Figure  46) 
reveals  that  the  strength  of  the  section  is  reduced  further  when 
subjected  to  pure  rotation  with  an  increase  in  the  r/t  ratio.  Also, 
an  interesting  point  to  note  is  that  the  strength  of  Section  CCS-2 
is  less  when  subjected  to  displacement/ rotation  ratio  5.0  compared 
to  an  imposed  pure  rotation  case,  while  the  strength  is  higher  for 
Imposed  displacement/rotation  ratio  of  15  compared  to  the  imposed 
rotation  case.  A reverse  phenomenon  Is  observed  in  the  case  of 
CCS-1  (Figure  46).  Figures  49.1  through  49.6  Indicate  the  M-0  and 
P-6  curves  for  various  deformation  ratios. 
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COMPRESSIVE  DISPLACEMENT  5 (Ins) 


Figure  47.1  CCS-1  Pure  Compression 
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Figure  47.2  CCS-1  Pure  Rotation 
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Figure  47.5  CCS- 1 Compresslon/Rotatlon  ■ 5.0 
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Figure  48  CCS  #2  - Failure  Envelope 
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Figure  49.1  CCS-2  Pure  Compression 
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F19MT*  «t.4  CCS- 2 £*U*t1< 


Figure  49.5  CCS- 2 Compress lon/Rotitlon  * 5.0 


(3)  Response  Of  CCS-3 


The  cylindrical  sections  analyzed  CCS-1,  CCS-2  seemed  to 
indicate  that  a change  in  the  r/t  ratio  changed  the  form  of  the 
failure  envelope.  Both  these  sections  were  subjected  to  imposed 
deformation  of  displacements  and  rotations.  The  response  of  a 
cylindrical  section  was  now  analyzed  when  subjected  to  a pressure 
load  of  ”3-4-5"  (Figure  50).  This  pressure  loading  was  idealized 
as  shown  in  Figure  50.  The  cylinder  was  symmetrically  loaded  with 
respect  to  vertical  and  horizontal  planes  through  the  longitudinal 
axis.  Such  a symmetric  loading  pattern  causes  symmetric  cracking 
and  yield  patterns.  Thus  only  a quarter  of  the  cylinder  was 
modeled  with  appropriate  boundary  conditions.  The  boundary 
conditions  imposed  on  the  section  were  as  follows:  symmetry  boundary 
condition  In  X-direction  on  face  (1-4-12-9),  symmetry  boundary 
condition  on  face  (9-12-24-21)  in  the  Y-direction  while  the  face 
(1-4-16-13)  was  restrained  in  the  Z-direction.  The  finite  element 
model  is  shown  in  Figure  51.  As  before,  concrete  was  modeled  using 
three-dimensional  brick  elements  while  steel  was  modeled  as  3-D 
truss  elements.  The  uniaxial  stress-strain  relations  for  concrete 
and  steel  are  shown  in  Figure  17  and  Figure  18.  The  ratio  r/t  was 
assumed  to  be  4. 

Figure  52  shows  a plot  of  load  vs.  displacement  at  the  crown 
of  the  cylinder  (node  4)  when  subjected  to  a monotonical ly  increasing 
distributed  load  of  "3-4-5".  As  can  be  seen  from  this  figure,  the 
section  yields  completely  in  isotropic  compression  at  point  2,  while 
failure  of  the  section  occurs  at  point  3 at  a displacement  of 
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Actual  cylindrical  section  with  actual  "3-4-5'  pressure 
load  distribution 


Actual  cylindrical  section  with  idealized  pressure  loading 

Figure  50.  Actual  and  Idealized  Loading.  Cylindrical  Section 
with  a Pressure  Load  of  "3-4-5" 


Figure  51.  Finite  Element  Model  for  Cylindrical  Section  with  a Pressure  Load  of  "3-4-5 
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SECTION  IV 

APPLICATIONS  OF  PROPOSED  MODEL 

1.  FAILURE  DESIGN  CURVES 

Based  on  the  parametric  data  of  Section  III,  failure  design  curves  are 
now  proposed  shown  in  Figures  53  through  56  for  rectangular  sections  while 
Figure  57  shows  the  curve  for  cylindrical  sections.  In  formulating  each  of 
these  sets  of  design  curves,  only  one  parameter  Is  varied  while  keeping  all 
the  other  parameters  constant.  As  an  example.  In  Figure  53,  only  the  uni- 
axial compressive  strength  of  concrete  is  varied  while  the  amount  of 
longitudinal  and  transverse  steel  Is  kept  constant.  It  should  be  emphasized 
that  the  dimensions  of  the  section  are  kept  constant. 


In  Figure  53  it  can  be  seen  that  the  maximum  deformation  of  the  section, 
at  failure,  reduces  as  the  strength  of  concrete  Increases.  Thus  a section  with 
concrete  having  an  ultimate  uniaxial  compressive  strength  of  6.0  ksi  has 
less  deformation  strength  than  the  one  with  a concrete  strength  of  5.0  ksi. 

In  both  cases,  the  ultimate  strain  of  the  uniaxial  compressive  stress-strain 
relation  was  0.003  in/in.  This  difference  in  strengths  could  be  due  to  the 
nature  of  uniaxial  stress-strain  relations  used  in  both  cases.  In  the  case 
of  concrete  with  an  ultimate  compressive  stress  of  6.0  ksi,  unloading  corres- 
ponding to  the  descending  portion  of  the  curve,  from  a strain  of  0.002  in/in 
to  a strain  of  0.003  in/in  is  at  a faster  rate  than  for  the  similar  portion 
of  the  curve  in  the  case  of  concrete  with  an  ultimate  uniaxial  stress  of 
5.0  ksi.  This  imples  that  failure  is  accelerated  in  case  of  the  higher 
strength  concrete  as  the  stress  state  corresponding  to  the  ultimate  uniaxial 
compressive  stress  is  reached  in  the  concrete  loaded  in  compression.  A 
refinement  of  the  failure  envelope  for  a higher  strength  concrete  may  be 
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obtained  with  a refinement  of  the  finite  element  mesh.  In  the  same 
Figure  53,  for  the  section  with  a uniaxial  ultimate  compressive  stress 
f'c  of  9.3  ksl,  the  maximum  deformation  at  failure  is  higher  than  the 
maximum  deformation  at  failure  of  sections  with  concrete  strengths  in 
uniaxial  compression  of  5.0  ksi  and  6.0  ksi.  This  difference  In  maximum 
deformations  is  exhibited  due  to  the  higher  ultimate  uniaxial  compressive 
strain  of  0.0045  in/in  used  in  the  case  of  sections  with  concrete  strength  in 
uniaxial  compression  of  9.3  ksi. 

Figure  54  shows  the  failure  design  curves  varying  f'c  of  the  section. 

The  longitudinal  reinforcement  is  kept  constant  (p  = p1  = 1/4%)  as  also 
the  transverse  reinforcement  (A$^  = 0.6  in2).  This  figure  indicates  that 
as  the  percentage  of  longitudinal  reinforcement  is  decreased,  the  section 
behaves  close  to  a plain  concrete  section.  This  figure  also  indicates  that 
given  the  same  ultimate  strain  for  concrete  having  an  ultimate  stress  of 
5.0  ksi  and  6.0  ksi,  the  section  with  higher  concrete  has  a higher  deformation 
strength. 

Design  failure  curves  for  sections  having  the  same  strength  concrete,  but 
variation  in  longitudinal  reinforcement,  are  shown  in  Figure  55.  Note  that 
the  ratio  of  longitudinal  to  transverse  reinforcement  was  constant  in  all  the 
three  cases.  It  can  be  seen  from  this  figure  that  the  section  with  1%  longi- 
tudinal reinforcement  has  more  deformation  strength  than  one  with  2% 
reinforcement  at  each  face.  It  should  be  mentioned  here  that  if  the  area 
of  transverse  reinforcement  would  have  been  kept  the  same  for  the  section  with 
1%  longitudinal  reinforcement  as  that  of  the  section  with  2%  reinforcement 
(Asj,  ■ 3.5  in2),  the  section  would  gain  additional  strength  and  ductility. 

This  phenomenon  has  been  demonstrated  in  Figure  56. 
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Figure  56  shows  the  failure  curves  of  a section  having  the  same  strength 


concrete,  percentage  of  longitudinal  reinforcement,  but  different  areas  of 
transverse  reinforcement.  These  curves  indicate  a definite  Increase  in  the 
strength  and  ductility  of  the  section  with  an  increase  In  transverse 
reinforcement. 

Oue  to  the  limited  study  performed  on  cylindrical  sections,  only  one 
set  of  design  curves  is  presented.  The  strength  of  concrete  was  the  same 
(5.0  ksi),  the  percentage  of  longitudinal  and  hoop  steel  was  constant  but 
the  r/t  ratio  (radius  to  thickness  ratio)  was  varied.  Such  a set  of  design 
curves  is  shown  in  Figure  57.  As  can  be  seen  from  this  figure,  as  r/t  is 
increased  from  4 to  8,  the  deformation  strength  of  the  section  increases  in 
the  compression  zone  while  it  decreases  in  the  tension  zone.  Near  the  pure 
rotation  zone,  a section  with  a smaller  r/t  ratio  has  a higher  value  of 
maximum  deformation  at  failure. 

It  should  be  mentioned  that  these  design  curves  are  preliminary  design 
curves  based  on  the  limited  parametric  study  undertaken.  A more  extensive 
study  is  required  to  verify  and  confirm  the  results  presented,  which  is 
beyond  the  scope  of  this  report. 

2.  AN  EXAMPLE  APPLICATION 

An  attempt  is  now  made  to  define  an  example  application  of  the  failure 
design  curves.  Different  modeling  techniques  to  obtain  deformation  quantities 
for  this  example  will  also  be  reviewed.  It  should  be  emphasized,  at  the 
outset,  that  while  obtaining  deformation  quantities  for  this  example  only 
nonlinear  analysis  techniques  are  to  be  used.  The  failure  design  curves 
cannot  be  used  if  the  deformation  quantities  have  been  obtained  using  a 
linear  analysis. 


I—-  " 
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Rectangular  Section.  Design  Failure  Curves 
Variation  in  p = p’.  Constant  f'  = 5.0  ksi 


Strain  (10'  in/in) 


Fig.  56  Rectangular  Section.  Design  Failure  Curves. 

Variation  of  Transverse  Reinforcement. 
Constant  f'  a 5.0  ksi,  p * p'  * 1/4t 
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Consider  a cylindrical  section  having  a specified  r/t  ratio  and  a 
specified  finite  length  of  the  section,  L,  as  shown  in  Figure  58.  Various 
deformation  quantities  for  this  section  can  be  obtained  by  modeling  the 
structure  appropriately.  If  one  were  interested  in  the  axial  deformations  of 
the  cylinder  when  subjected  to  forces  or  deformations  axially,  one  could  model 
the  structure  as  a series  of  discrete  springs  with  lumped  masses  at  appropriate 
intervals  (Figure  59).  The  system  could  be  modeled  to  have  only  axial  degrees 
of  freedom;  i.e.,  the  masses  have  no  rotational  inertia  associated  with  them. 
Depending  on  the  type  of  problem  to  be  solved,  appropriate  boundary  conditions 
could  be  imposed  on  this  model.  The  model  is  then  analyzed  for  the  axial 
imposed  forces  or  deformations  and  the  distribution  of  axial  deformation  due 
to  these  imposed  forces  calculated.  Once  the  axial  deformations  are  known 
one  could  go  back  to  the  design  failure  curves  and  check  the  deformations 
(converted  to  appropriate  strain  quantities)  against  the  axial  deformations 
the  section  can  resist  as  indicated  by  the  strain  axis  of  the  failure  envelope. 

If,  however,  the  maximum  rotation  attainable  at  any  section  of  the  cylinder 
was  to  be  checked  at  various  points  along  the  length  of  the  section  when  Imposed 
with  certain  moments  or  rotations,  with  no  axial  forces,  one  could  model  the 
same  structure  as  a series  of  beams  interconnected  at  various  Intervals  along 
the  length  of  the  section  (Figure  60). 

Concentrated  plastic  hinges  could  be  formed  at  the  end  of  these  sections. 
Rotation,  as  well  as  extension  or  compression  of  these  plastic  hinges,  could 
be  permitted.  If  the  plastic  hinges  had  only  rotation  associated  with  them; 
then  one  would  have  a case  of  pure  rotation  at  the  displacement  degrees  of 
freedom.  One  could  check  these  maximum  rotations  so  obtained  against  the 
rotational  strength  given  by  the  failure  envelope  of  the  section. 
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Axial  imposed 
forces  or  imposed 
deformations 


Only  axial  degrees  of  freedom  considered 


* 

Cylindrical  Section  Modeled  as  Discrete  Springs 
and  Lumped  Masses 


Consider  the  cylinder  being  subjected  to  a pressure  loading  of  3-4-5  as 


indicated  in  Figure  61.  Local  buckling  of  such  a section  could  be  an 
important  mechanism  of  failure.  One  could  model  such  a structure  as  shown 
in  Figure  61.  In  this  figure,  it  can  be  seen  that  the  transverse  section  of 
the  cylinder  is  modeled  with  beam  elements.  A critical  portion  of  the  cylinder 
can  be  isolated  (Figure  61)  and  can  be  modeled  as  beam  elements  connected 
together  by  hinges.  Appropriate  boundary  conditions  can  be  applied  at  the 
hinges  and  the  deformations  obtained  from  a nonlinear  analysis.  Hinge 
behavior  is  graphically  shown  in  Figure  61,  where  these  hinges  are  capable 
of  both  elastic  and  plastic  displacements  and  rotation.  Note  that  in  this 
case,  the  deformation  capacity  would  be  checked  against  the  deformation  of 
the  design  failure  envelopes  for  the  rectangular  section,  rather  than  the 
cylindrical  section. 

All  the  above  cases  of  modeling  of  structures  and  their  application  with 
respect  to  design  failure  curves  fonn  the  analytical  portion  of  the  application 
of  design  failure  curves.  One  could  also  use  experimental  results  to  verify 
these  design  failure  curves. 

3.  FAILURE  ENVELOPES  AND  YIELD  ENVELOPES 

The  failure  envelopes  as  obtained  from  a three-dimensional  finite  element 
analysis,  were  presented  in  dimensionless  quantities  of  strain  and  curvature 
in  Section  III  of  this  report.  These  failure  envelopes  were  reduced  to 
design  failure  envelopes  in  terms  of  strains  and  curvatures.  It  should  be 
emphasized  here  that  the  strains  and  curvatures  presented  In  these  curves 
are  referenced  at  the  geometric  centroid  of  the  section.  Thus  the  strain 
at  failure  is  the  extensional  or  compression  strain  at  the  geometric  centroid 
of  the  section  while  the  curvature  at  failure  is  the  curvature  at  the 


and  plastic  displacements  and 
rotations  when  subjected  to 
pressure  loading. 


P 


Hinge  Behavior 


Figure  61.  Cylindrical  Section  Modelled  as  Beam  Elements. 


Strain  (10  in/in) 


Compression 


Fig.  62  CRS  #3.  Failure  Envelope,  Ultimate, 
Envelope  and  Yield  Envelope 


geometric  centroid  of  the  section.  Note  that  at  failure  the  strains  in  the 
extreme  fibers  are  different  from  the  strains  presented  In  the  failure 
envelopes  and  in  the  design  failure  envelopes  of  each  section.  To  measure 
the  ductility  of  the  section,  an  Initial  yield  surface  has  to  be  defined 
compared  to  the  failure  envelope  of  the  section.  The  yield  surface  was 
defined  when  the  effective  stress  at  any  point  in  the  section  exceeded  a 
value  of  tq  which  was  a function  of  fc»  *c.  ft»  ^bc  where  f^.  indicates  the 
ultimate  strength  In  uniaxial  compression,  while  fc,  ft,  and  f^  indicate 
initial  yield  strength  in  uniaxial  compression,  tension  and  equal  biaxial 
loading,  respectively.  The  initial  yield  strength  in  uniaxial  compression, 
fc,  was  assumed  to  be  0.65fc-  A yield  surface  defined  in  such  a way  is 
shown  for  a single  section  CRS-3  in  Figure  62. 


4.  FAILURE  ENVELOPE  AND  ULTIMATE  ENVELOPE 

A comparison  between  the  failure  envelope  obtained  from  a 3-D  finite 
element  analysis  with  a nonlinear  constitutive  material  model  and  the  ultimate 
strength  envelope  using  a one  dimensional  fiber  analysis  is  shown  in  Figure  62. 
The  comparison  of  these  two  types  of  envelopes  refers  to  CRS  #3. 

It  is  appropriate  at  this  stage  to  reiterate  the  difference  in  loading 
of  the  section  to  obtain  failure  strength  using  a 3-D  FEA  analysis  as  compared 
to  the  ultimate  strength  as  obtained  from  a one  dimensional  fiber  analysis. 

In  the  analysis  of  a failure  strength  of  a section,  using  a 3-D  FEA,  the 
neutral  axis  of  imposed  deformations  is  fixed  while  the  deformations  are 
monotonically  Increased  in  a specific  ratio  of  axial  to  rotational,  until  the 
failure  strength  of  the  section  is  reached.  The  section  was  subjected  to 
seven  specific  fixed  deformation  ratios  to  obtain  seven  failure  strengths  to 
give  a failure  envelope  for  that  section.  In  contrast,  in  a fiber  analysis. 
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the  neutral  axis  changes  position  within  the  section,  as  the  imposed 
strain  distributions  are  varied  from  zero  to  a distribution  at  which  the 
ultimate  strength  of  the  section  is  reached. 


The  similarity  between  the  failure  envelope  and  the  ultimate  strength 
envelope  (P-M  curve.  Figure  63)  should  also  be  emphasized.  Each  point  on 
the  failure  envelope  has  a unique  position  of  neutral  axis  of  Imposed 
deformation  associated  with  it.  Thus  all  points  In  the  failure  envelope 
have  a different  position  of  neutral  axis  of  Imposed  deformations.  A 
similar  phenomenon  is  observed  in  the  case  of  the  ultimate  strength  envelope 
obtained  from  a one  dimensional  fiber  analysis. 

Thus,  if  we  have  an  ultimate  envelope  defined  for  an  ultimate  strain 
of  0.002  In/in  in  the  extreme  fiber  of  concrete  in  compression,  then  with 
each  point  on  this  curve  is  associated  a certain  strain  distribution  with 
a location  of  the  neutral  axis  unique  to  that  point.  Thus,  If  we  consider 
a point  A on  this  ultimate  envelope  of  the  P-M  interaction  diagram,  then 
the  strain  distribution  associated  with  this  point  is  shown  in  Figure  63. 

From  the  strain  distribution,  the  curvature  and  hence  the  displacement  (6) 
and  rotation  (6)  for  the  section  can  be  easily  calculated.  Given  this  6 and 
9,  the  ratio  (6/0)  of  deformations  causing  the  ultimate  condition  can  be 
calculated.  This  6/9  ratio  compared  to  a similar  5/9  ratio  for  the  failure 
envelope  using  a three-dimensional  analysis  Is  shown  In  Figure  62.  Note 
that  a sample  calculation  for  such  a point  on  the  ultimate  envelope  Is 
shown  In  Appendix  B.  The  displacement  and  rotation  at  the  geometric  centroid 
of  the  section  Is  calculated  corresponding  to  the  balanced  point  on  the  P-M 
Interaction  diagram.  The  concrete  cover  was  considered  In  the  calculations 
which  was  neglected  in  the  development  of  the  failure  envelopes. 
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5.  LIMITATIONS  OF  THE  PROPOSED  MODEL 

While  developing  the  failure  envelopes  for  various  sections,  certain 
assumptions  were  made  in  the  analysis  which  will  be  summarized  in  this 
section. 

(1)  The  failure  envelopes  were  presented  in  the  terms  of  strains  and 
curvatures  for  both  the  rectangular  and  the  cylindrical  sections.  It  should 
be  emphasized  that  these  strains  and  curvatures  were  referenced  to  the 
geometric  centroid  of  the  section.  Thus,  the  strain  was  the  extension  or 
compressive  strain  at  the  geometric  centroid  of  the  section  while  the 
curvature  was  the  curvature  at  the  centroid  of  the  section. 

(2)  The  design  failure  curves  were  presented  in  three  different  forms 
for  the  rectangular  section:  One,  failure  curves  obtained  with  variation  of 
the  ultimate  uniaxial  compressive  stress  of  concrete  f£;  two,  failure  curves 
obtained  by  varying  the  percentage  of  longitudinal  reinforcement;  and  three, 
curves  obtained  by  varying  the  transverse  reinforcement.  By  doing  so,  upper 
and  lower  bounds  were  obtained  for  each  of  the  set  of  curves.  Thus,  for  the 
first  set  of  curves,  the  upper  bound  was  defined  by  an  ultimate  uniaxial 
compressive  stress  of  9.0  ksi,  while  the  lower  bound  was  defined  by  concrete 
having  an  ultimate  uniaxial  compressive  stress  of  5.0  ksi.  Note  that  the 
failure  curves  are  very  sensitive  to  the  type  of  uniaxial  compressive  stress 
strain  curves  for  concrete  Input  to  the  model.  The  upper  and  lower  bound 
of  percentage  of  reinforcement  is  defined  as  2%  and  1/4%  of  reinforcement 
at  each  face.  Note  also  the  sensitivity  of  the  curves  to  the  variation  of 
this  parameter.  The  section  with  1/4%  reinforcement  on  each  face  behaves 
almost  like  a plain  concrete  section. 

(3)  Only  two  basic  types  of  geometric  sections  were  considered:  the 


rectangular  cross-section  and  the  cylindrical  cross-section.  The  dimensions 
of  these  sections  were  kept  constant  for  all  the  analyses.  Note  that  the 
effect  of  variation  in  dimensions  on  the  failure  envelope  was  not  studied. 

Also,  the  concrete  cover  for  the  rectangular  and  the  cylindrical  section  was 
neglected  in  the  analysis.  It  is  believed  that  such  an  assumption  would  not 
influence  the  failure  curves  significantly. 

(4)  Deformations  instead  of  forces  were  Imposed  on  the  section.  These 
deformations  were  imposed  as  a combination  of  displacements  and  rotations  in 
a certain  ratio  which  was  kept  fixed  for  a particular  analysis.  Thus,  for  a 
given  analysis  that  gave  a point  on  the  failure  envelope,  the  neutral  axis  of 
deformation  was  fixed  while  the  deformations  were  monotonical ly  increased  until 
failure  of  the  section  occurred.  This  assumption,  it  is  believed,  does  not 
affect  the  final  failure  deformations  obtained  from  analyses.  However,  care 
should  be  taken  to  Interpret  the  results  obtained  up  to  failure  deformations; 
i.e.,  the  M-e  and  P-5  curves  obtained  from  such  an  analysis. 

(5)  It  should  be  noted  that  only  flexural  and  axial  deformations  were 
considered  in  the  development  of  these  failure  design  envelopes.  The  effect 
of  shear  deformations  on  these  failure  envelopes  was  neglected.  The  effect 
of  shear  and  the  interaction  failure  envelopes  of  shear,  axial  and  flexural 
deformations  form  a different  aspect  of  the  problem  which  has  to  be  considered 
as  a separate  item  of  research. 

(6)  Certain  engineering  judgments  will  have  to  be  made  while  using  these 
failure  design  curves.  Note  that  only  symmetrically  reinforced  sections  were 
considered  in  the  analysis.  Thus,  if,  for  example,  one  were  considering  a 
rectangular  section  which  was  unsymmetrically  reinforced,  say  1*  on  the  top 
face  and  2%  at  the  bottom  face,  then  one  would  expect  the  failure  envelope  to 
be  different  in  the  compression  and  the  tension  zone.  The  curve  may  not  be 
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very  different  in  the  tension  zone,  where  the  failure  deformation  for  an 
imposed  deformation  of  extension/rotation  = 5.0  is  considered.  This  type  of 
behavior  would  be  expected  since  the  compressive  zone  of  the  section  is  small, 
and  steel  In  the  tension  zone  of  the  section  controls  behavior.  If,  on  the 
other  hand,  a similar  unsymmetrically  reinforced  section  was  considered  to 
establish  failure  deformations  for  an  imposed  deformation  pattern  of 
compression/rotation  of  5.0,  a slight  increase  in  the  section  capacity  may 
be  expected,  compared  to  a section  reinforced  symmetrically  with  2%  rein- 
forcement on  each  face.  This  would  be  due  to  the  fact  that  the  ductility  of 
the  section  can  be  utilized  to  a larger  extent.  In  this  case,  where  the 
section  is  reinforced  symmetrically  with  2%  reinforcement  on  each  face, 
before  the  yield  of  re-bars  in  compression,  it  Is  this  steel  in  compression 
which  resists  a large  portion  of  the  load.  As  soon  as  this  steel  yields, 
the  redistribution  of  stresses  in  the  concrete  causes  failure  of  the  concrete 
section  at  a smaller  displacement  than  if  the  section  was  symmetrically  reinforced 
with  1%  reinforcement  in  each  face.  Thus,  with  an  unsymmetrically  reinforced 
section  with  1%  steel  on  the  top  face  and  2%  steel  on  the  bottom  face,  the 
failure  deformation  for  an  imposed  deformation  of  compression/rotation  of 
5.0  would  be  somewhere  in  between  the  failure  deformation  obtained  if  the 
section  was  symmetrically  reinforced  on  each  face  with  1%  and  2%  longitudinal 
reinforcement. 

(7)  The  failure  strengths  and  failure  envelopes  are  very  sensitive  to 
type  of  uniaxial  compressive  stress-strain  curve  for  concrete.  Thus,  In  using 
the  failure  design  curves,  to  compare  maximum  deformations  of  critical  sections 
of  a structural  system  obtained  using  a nonlinear  analysis  (technique  outlined 
in  Section  IV  as  an  example  application  of  the  design  curves),  it  is  essential 
that  material  behavior  of  concrete  in  uniaxial  compression  of  this  structural 


system  be  similar  to  the  o-e  curves  used  to  obtain  failure  envelopes  for 
various  sections. 


(8)  In  conclusion,  it  must  be  emphasized  that  while  using  these  failure 
design  curves,  extreme  caution  should  be  taken.  The  assumptions  involved  in 
their  development  should  be  thoroughly  studied.  Any  given  problem  that  is 
being  analyzed  should  be  reduced  to  a problem  similar  to  the  ones  studied  in 
this  report  before  these  curves  are  applied.  Certain  engineering  judgments 
may  have  to  be  made  in  the  process. 


APPENDIX  A 


1.  INTRODUCTION 

A brief  discussion  and  explanation  of  the  program  is  presented  in 
this  section.  The  organization  of  the  auxiliary  subprogram,  the  three- 
dimensional  8 to  20  node  finite  element  solution,  and  the  material  con- 
stitutive model  is  discussed  herein.  For  a complete  summary  of  the  base 
program,  the  reader  Is  referred  to  Reference  4.  A complete  user's  guide 
is  also  presented  at  the  end  of  this  appendix. 

2.  PROGRAM  LAYOUT  AND  CAPABILITIES 

The  main  program  ANSR  is  set  up  in  two  basic  parts;  the  base  program 
and  the  auxiliary  programs.  Each  auxiliary  program  Is  a set  of  subroutines 
defining  a specific  element  type.  Each  auxiliary  program  has  at  least  four 
main  subroutines. 

1.  INEL:  Input  and  initialization  of  element  information. 

2.  STIF:  Formation  of  element  tangent  stiffness  in  static  analysis, 

or  effective  element  stiffness  in  dynamic  analysis. 

3.  RESP:  Computation  of  response  of  the  elements  to  Imposed  deforma- 

tions or  imposed  forces. 

4.  OUT:  Output  of  envelope  values  of  element  deformation  and  actions 

at  specified  load  increments  in  static  analysis  or  at 

specified  time  intervals  in  dynamic  analysis. 

Each  of  these  routines  is  identified  by  a number.  Two  types  of  elements 
exist  in  the  program,  the  three-dimensional  truss  elements  identified  as 
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element  number  1,  and  the  three-dimensional  8 to  20  node  finite  element, 
specified  as  element  number  4.  Hence,  subroutines  INEL1,  STIF1,  RESP1 , 

0UT1  all  refer  to  elanent  number  1,  while  subroutines  INEL4,  STIF4,  RESP4, 
0UT4  refer  to  element  number  4.  Each  of  these  main  subroutines  could  refer 
to  numerous  secondary  subroutines  defining  material  behavior. 

The  elastic  plastic  fracture  strain  hardening  material  model  as  proposed 
by  Chen  and  Chen  with  modifications  was  identified  as  material  Model  5.  The 
subroutines  describing  the  constitutive  relations  are  discussed  herewith 
with  an  explanation  of  some  of  the  parameters  and  the  significance  of  the 
values  assigned  to  some  of  the  parameters.  There  are  12  secondary  sub- 
routines that  define  the  material  behavior.  A brief  description  of  each  of 
these  subroutines  Is  presented. 

1.  ADRS45:  This  subroutine  assigns  addresses  to  the  material  data 

In  the  element  information  array. 

2.  CMAT45:  This  subroutine  evaluates  the  elasticity  matrix  for  the 

material  model.  For  details  of  the  matrix  refer  to 
Section  II. 2 of  this  report. 

3.  EPSLOP:  This  subroutine  evaluates  the  strain  hardening  function  H. 

4.  CMAT45A: -This  subroutine  forms  the  elasticity  matrix  for  the 

material  model. 

5.  YVALUE:  This  subroutine  evaluates  the  value  of  effective  stress 

using  the  appropriate  loading  function,  after  the  stress 
state  has  been  established  as  being  either  a compressive 
stress  state,  tension-compression  stress  state  or  a 
tensile  stress  state.  For  details  of  the  loading  func- 
tion, refer  to  Section  II. 2.  The  loading  function  is 
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found  to  be  inconsistent  with  the  criteria  to  check 
the  stress  state  for  certain  conditions  of  loading 
(Section  II. 2).  Thus,  If  the  stress-state  established 
is  a pure  compressive  one  and  the  effective  stress  so 
calculated  found  to  be  zero,  the  stress  state  so 
established  Is  inconsistent.  Hence,  it  Is  now  assumed 
to  be  a tension  compression  stress  state.  All  the 
requisite  parameters  are  reset  and  the  effective  stress 
calculated  using  the  appropriate  loading  function.  Note 
that  the  criteria  to  establish  the  stress-states  are 
inconsistent  only  in  the  compression  and  tension-compression 
region. 

TWXIJ:  This  subroutine  evaluates  the  first  invariant  of  stress  I, 

and  second  invariant  of  devlatoric  stress  J,. 

0SHT45:  This  subroutine  determines  the  overshoot  Into  the  plastic 
state.  The  overshoot  is  determined  by  the  roots  of  the 
quadratic  equation  (Section  II. 2,  Equation  18).  Both  roots 
are  evaluated  and  the  minimum  positive  root  Is  the  required 
value. 

MANV45:  This  subroutine  prints  the  peak  stress  components  and  strain 
components  for  each  element. 

PRIN45:  This  subroutine  prints  the  values  of  stress  components 
and  strain  components  for  each  element  at  the  end  of 
specified  load  steps  for  static  analysis  and  at  specified 
time  Intervals  for  dynamic  analysis. 
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10.  MATP45:  This  subroutine  reads  In  appropriate  material  data, 

evaluates  all  the  material  properties  and  parameters, 
and  prints  the  material  property  table. 

11.  MRES45:  This  subroutine  is  the  routine  which  evaluates  the 

material  response.  State  determination  phase  occurs 
at  this  level  in  this  subroutine.  Appropriate  stress 
states  are  established  following  which  the  crack  behavior 
or  the  yield  behavior  Is  evaluated.  The  stress  incre- 
ment is  calculated  assuming  elastic  behavior  for  a given 
Increment  of  loading  or  deformation.  The  stress  Is  then 
checked  for  being  outside  the  initial  discontinuous  surface. 
If  it  has  overshot  into  the  plastic  range,  then  the  elastic 
increment  is  appropriately  scaled  down  and  the  remaining 
portion  of  the  plastic  strains  evaluated  (Section  II. 2, 
Equation  22).  The  codes  for  each  type  of  behavior  are  set 
in  this  subroutine. 

12.  STRAIN  This  subroutine  evaluates  the  equivalent  plastic  strain. 

This  subroutine  is  called  at  various  stages  by  MRES45. 

The  labelled  common  blocks  used  in  the  above  12  subroutines  are: 

1.  C0MM0N/C0NCRT/AAU(2,10) ,AA0(2,10) ,TAU(2,10) ,TA0( 2,10) ,AL FA( 2,10) , 

BETA ( 2,10) ,XK(2,10),XN(2,10) ,PHC (10,10) , PHT (10,10) , 
PFC(10, 10), PFT (10,10) ,STRN0(10) , STRNF (10) ,C0EFF. 
where  AAU,  AA0  = Ay  and  Aq  of  Section  1 1. 2 

Equations  5,  6,  9,  10. 

TAU,  TA0  = tu  and  tq  of  Section  II. 2 
Equations  7,  8,  11 , 12. 
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ALFA,  BETA  = a,  8 of  Section  II. 2 
Equations  3,  4. 

XK,  XN  = The  first  and  the  section  coefficient 

in  the  loading  functions  In  Section  II. 2, 
Equations  1,  2. 

PHC.PFC  = Strain  and  stress  for  the  nonelastic 

portion  of  the  equivalent  stress  strain 
curve  in  compression  region. 

PHT.PFT  = Strain  and  stress  for  the  nonelastic 

portion  of  the  equivalent  stress-strain 
curve  in  the  tension  region. 

STRNU  = Equivalent  ultimate  plastic  strain. 

STRNF  = Equivalent  fracture  plastic  strain 

COEFF  = The  coefficient  by  which  the  elasticity 
matrix  in  subroutine  CMAT45  is  multiplied 
once  the  gauss  integration  point  In  the 
element  is  either  cracked  or  crushed. 

Note  that  a maximum  of  10  different  material  types  can  be 

handled  by  the  program.  Also,  the  parameters  A^,  Aq,  t , tq,  a,  8, 

XK,  XN  are  different  in  the  compression  region  and  the  tension 
compression  region. 


COMMON/CR ITER/ CONSTA , CONSTB , YLDST , ULTST , GK , GN , I DTC , QFC ( 1 0 ) , KKCC . 
where  CONSTA  = a 

CONSTB  = B 

YLDST  = tq 

ULTST  = x 
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GK 


XK 


GN  = XN 

IDTC  = Indicator  for  establishing  stress 

states.  If  IDTC  = -1,  then  the 
stress  state  is  a tension-compression 
stress  state;  If  IDTC  = 1,  then  the 
stress  state  is  pure  compression.  If 
IDTC  = -6,  then  the  stress  state  is  a 
pure  tension  one. 

KKCC  = Number  of  nonelastic  portions  of 

the  uniaxial  stress-strain  curve 
(Figure  A.l ) . 

The  possible  codes  printed  at  the  end  of  specified  load  steps 
for  each  element  of  element  number  4,  for  static  analysis,  or  at 
specified  dynamic  analysis,  are  indicated  in  Figures  3,  4 and  5. 


i 
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A. 


PROBLEM  INITIATION  AND  TITLE  (A5,  3X,  18A4)  - One  card 
Columns  1 - 5:  Punch  the  word  START 
6 - 8:  Blank 

9 - 80:  Problem  title,  to  be  printed  with  output. 
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B.  NODE  INFORMATION 


Bl.  CONTROL  INFORMATION  (915)  - One  card 

Columns  1 - 5:  Total  number  of  nodes. 


6-10:  Number  of  "control"  nodes,  for  which 

coordinates  are  specified  directly 
(NCNOD).  See  section  B2. 

11  - 15:  Number  of  coordinate  generation  commands 
(NODGC).  See  section  B3. 

16  - 20:  Number  of  commands  specify inn  nodes 

with  zero  displacements  (NDCON).  See 
section  B4. 


21  - 25:  Number  of  commands  specifying  nodes  with 
egual  displacements  (NIDDOF).  See 
section  B5. 


26  - 30:  Number  of  commands  specifying  nodal 
masses  (NMSGC).  See  section  B6. 

31  - 35:  Number  of  element  orouDS  (NELGR,  max.  20). 
See  section  E. 

40:  Execution  code  (KEXEC)  as  follows: 

(a)  zero  or  blank:  full  execution. 

(b)  1:  data  checkino  only. 

( c )  - 1 : full  execution,  but  only  if  the 

structure  stiffness  and  element 
data  can  be  held  in  core. 

45:  Stiffness  storage  code  (KSCHM) , as  follows. 

(a)  zero  or  blank:  duplicate  stiffness 

matrix  held  in  core. 

(b)  1:  duplicate  stiffness  matrix  stored 

on  scratch  file. 

50:  Number  of  imposed  displacement  directions 

(NC0RD).  Maximum  =3.  If  specified  as  a 
negative  quantity,  then  a printout  of  axial 
force  P,  Moment  M and  unbalanced  load  is 
obtained  at  the  end  of  each  load  step. 
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B2.  CONTROL  NODE  COORDINATES  (I5.3F10.0)  - NCNOD  cards 


Columns  1 - 5:  Node  number.  In  any  sequence. 

6-15:  X coordinate. 

16-25  Y coordinate. 

26  - 35:  Z coordinate. 
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B3.  COORDINATE  GENERATION  (415, FIO. 0,1015)  - NODGC  cards. 


Note: 


Columns  1 - 5:  Node  number  at  beginning  of  generation 

line.  This  must  either  be  a control  node, 
or  must  have  been  generated  by  a previous 
generation  command. 


6-10:  Node  number  at  end  of  generation  line. 

This  node  must  also  have  been  specified 
previously. 

11  - 15:  Number  of  nodes  to  be  generated  along  line. 
If  the  nodes  to  be  generated  are  listed 
in  Columns  31-80,  this  number  may  not 
exceed  10. 

16  - 20:  Node  number  difference  between  successive 

generated  nodes,  and  between  first  generated 
node  and  node  at  beglnnlnq  of  generation 
line.  May  be  negative.  Leave  blank  If 
generated  nodes  are  listed  In  Columns  31-80. 


21  - 30:  Spacing  between  nodes,  as  follows: 

(a)  zero  or  blank:  generated  nodes  are 
spaced  uniformly  along  the  neneratlon 
line. 

(b)  less  than  1.0:  spacing  between  nodes 
is  this  proportion  of  the  length  of  the 
oeneration  line. 

(c)  1.0  or  larger:  spacing  between  nodes 
Is  equal  to  this  distance. 


31  - 80:  Up  to  10  fields,  each  15.  List  nodes  to 
be  generated.  In  sequence  along  generation 
line.  Required  only  If  Columns  1 6-20  are 
blank. 


I 


It  Is  not  necessary  to  provide  coordinate  generation  commands  for 
nodes  which  are  sequentially  numbered  between  the  beginning  and  end 
nodes  of  any  straight  line,  and  which  are  equally  spaced  along  that 
line.  After  all  generation  coitmands  have  been  executed,  the 
coordinates  for  each  group  of  unspecified  nodes  are  automatically 
generated  assuming  sequential  numbering  and  equal  spaclno  along 
a line  joining  the  specified  nodes  Immediately  preceding  and 
following  the  group.  That  Is,  any  generation  command  with  a node 
number  difference  of  one  and  equal  spacing  Is  superfluous. 
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B4. 


NODES  WITH  ZERO  DISPLACEMENTS  ( I5.4X.6I1 ,1 315)  - NDCON  cards 

Columns  1 - 5:  Node  number,  or  number  of  first  node  In  a 

series  of  nodes  covered  by  this  command.  See 
Note  following  for  repetition  of  nodes. 


10:  Constraint  code  for  X displacement,  as  follows. 

(a)  zero  or  blank:  displacement,  not  con- 
strained to  be  zero. 

(b)  1:  displacement  constrained  to  be  zero. 

11:  Code  for  Y displacement. 

12:  Code  for  Z displacement. 

13:  Code  for  XX  rotation. 

14:  Code  for  YY  rotation. 

15:  Code  for  ZZ  rotation. 

16  - 20:  Number  of  last  node  In  series  of  nodes 
covered  by  this  command.  Leave  blank 
or  punch  zero  for  a single  code,  or  If 
the  nodes  In  the  series  are  listed  In 
columns  31 -HO. 

21  - 25:  Node  number  difference  between  successive 
nodes  In  series.  Leave  blank  for  a single 
node,  or  If  the  nodes  In  the  series  are 
listed  In  Columns  31-80. 

29  - 30:  Number  of  nodes  listed  In  Columns  31-80, 
following.  This  list  Is  considered  only 
if  Columns  16-20  are  blank  or  zero.  Leave 
blank  for  a single  node. 

31  - 80:  Up  to  10  fields,  each  15.  List  second, 
etc.  node  of  series. 


Note:  If  constraint  codes  are  specified  more  than  once  for  any  node, 
the  last  specified  value  is  assumed.  For  plane  or  axlsymmetric 
problems,  the  first  command  should  cover  all  nodes  and  should 
constrain  all  except  the  relevant  displacements.  Additional 
cards  to  modify  the  constraint  codes  at  particular  nodes  should 
then  be  added. 
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B5. 


NODES  WITH  EQUAL  DISPLACEMENTS  (611 ,4X, 1415)  - NIDDOF  cards 


Columns 


1: 


2 

3 

4 

5 

6 


7 - 10 


Equal  displacement  code  for  X displacement, 
as  follows. 

(a)  zero  or  blank:  displacement  not  constrained 
to  be  identical. 

(b)  1:  displacement  for  all  nodes  in  group. 

Code  for  Y displacement. 

Code  for  Z displacement. 

Code  for  XX  rotation. 

Code  for  YY  rotation. 

Code  for  ZZ  rotation. 

Blank 


11  - 15:  Number  of  nodes  in  group. 

16  - 80:  Up  to  13  fields,  each  15.  List  nodes  in 
group.  The  first  node  must  be  the  small- 
est numbered  node  In  the  aroup.  See  Note 
following. 


Note:  If  the  group  has  more  than  thirteen  nodes,  specify  the  remaining 
nodes  on  additional  equal  displacement  commands.  The  smallest 
numbered  node  In  the  group  must  be  the  first  node  In  the  list 
for  all  comnands  defining  the  group.  Greater  computational 
efficiency  may  be  obtained  by  constraining  nodes  to  have  equal 
displacements.  However,  the  effect  of  specifying  equal  displace- 
ments may  be  to  Increase  the  band  width  of  the  structure  stiff- 
ness matrix.  This  may  result  in  an  Increase  in  the  required 
stiffness  matrix  storage  and/or  the  computational  effort  required 
to  solve  the  equations  of  motion.  Equal  displacements  specifications 
should  therefore  be  used  with  caution.  It  should  be  noted  that  the 
equation  solver  used  in  the  program  Is  less  sensitive  to  local 
Increases  in  the  stiffness  matrix  band  width  than  a conventional 
equation  solver  based  on  a banded  storage  scheme. 
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B5.1  IMPOSED  DISPLACEMENT  SPECIFICATION  - NCORD  sets  of  cards.  (See 
Section  bl ). 

Each  set  consists  of  an  imposed  displacement  direction  card  and  a 
node  specification  card  specifying  the  nodes  at  which  the  imposed 
displacements  are  applied. 

B5.1 (a)  IMPOSED  DISPLACEMENT  DIRECTION  SPECIFICATION  (4I5,?F10.0)  - One  card. 
Columns  1 - 5:  Specified  displacement  directions. 


(a) 

1:  X-direction. 

(b) 

2:  Y-direction. 

(c) 

3:  Z-direction. 

Specified  lever  arm  axis. 

(a) 

1:  X-axis. 

(b) 

2:  Y-axis. 

(c) 

3:  Z-axis. 

11  - 15: 

Number  of  cards  to  describe  nodes  defining 
neutral  axis  (NDC). 

16  - 20: 

Number  of  cards  to  describe  other  nodes 
(rotation  + displacement  nodes)  off  neutral 
axis  (NRC). 

21  - 30: 

Specified  displacement  Increment. 

31  - 40: 

Specified  rotation  (of  neutral  axis)  increment. 

85.1(b)  NODE  DESCRIPTION  (315)  - NDC  cards  followed  by  NRC  cards. 

(See  Section  B5.1 (a)) 

Columns 

1 - 5: 

First  node  in  the  series. 

6 - 10: 

Last  node  in  the  series. 

11  - 15: 

Node  increment. 

B6. 


NODAL  MASSES  ( 15, 6F10. 0,215)  - NMSGC  cards 


Columns  1 - 5: 

6 - 15: 
16  - 25: 


Node  number,  or  number  of  first  node 
in  a series  of  nodes  covered  by  this 
command. 

Mass  associated  with  X-dlsplacement 
degree  of  freedom. 

Mass  associated  with  Y-displacement 
degree  of  freedom. 


26  - 35:  Mass  associated  with  Z-dlsplacement 
degree  of  freedom. 

36  - 45:  Mass  associated  with  X-rotation 
deqree  of  freedom. 

46  - 55:  Mass  associated  with  Y-rotation 
degree  of  freedom. 

56  - 65:  Mass  associated  with  Z-rotation 
degree  of  freedom. 

66  - 70:  Number  of  last  node  In  series  of  nodes 
covered  in  series.  Leave  blank  for 
a single  node. 


Note: 


The  specification  commands  for  lumped  masses  will  generally  permit 
the  user  to  input  the  nodal  masses  with  only  a few  data  cards. 

Any  node  may,  if  desired,  appear  in  more  than  one  specification 
command.  In  such  cases  the  mass  associated  with  any  degree  of 
freedom  will  be  the  sum  of  the  masses  specified  in  separate 
commands.  If  certain  nodes  are  constrained  to  have  an  equal 
displacement,  the  mass  associated  with  this  displacement  will  be 
the  sum  of  the  masses  specified  for  the  Individual  nodes.  If  a 
mass  is  specified  for  any  degree  of  freedom  that  is  constrained 
to  be  zero,  it  is  Ignored. 
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C.  LOAD  SPECIFICATION 

Cl.  CONTROL  CARD  (10I5.3F10.0)  - One  card 


Columns  1 - 5:  Code  for  static  and/or  dynamic  analysis, 

(KSTAT) . 

(a)  1:  static  analysis  only. 

(b)  2:  thermal  analysis  only. 

(c)  3:  dynamic  analysis  only. 

(d)  6:  static  analysis  followed  by 

dynamic  analysis. 

6-10:  Number  of  static  force  patterns  to  be 

specified  (NSPAT).  See  section  C2. 

If  blank  or  zero,  no  static  loads  will 
be  applied. 

11  - 15:  Number  of  static  force  application  commands 
(NSLGC).  See  Section  F. 

16  - 20:  Number  of  thermal  load  patterns  to  be 
specified  (NTPAT). 

21  - 25:  Number  of  thermal  load  application  commands. 

26  - 30:  Code  for  ground  motion  records  (IGM),  as 

follows: 

(a)  zero  or  blank:  no  ground  motion  records. 

(b)  1:  ground  motion  records  will  be 

specified.  See  Section  C3. 

31  - 35:  Number  of  dynamic  force  records  to  be 
specified  (NDLR).  See  Section  C4. 

36  - 40:  Largest  number  of  points  on  any  dynamic  force 
record.  This  number  is  used  for  storage 
allocation. 

41  - 45:  Number  of  commands  defining  points  of  application 
of  dynamic  force  records  (NDLGC).  See  Section  C5. 

46  - 50:  Number  of  integration  time  steps  to  be  considered 
in  dynamic  analysis. 

51  - 60:  Integration  time  step.  At. 

61  - 70:  Integration  method  parameter,  6,  in  Newmark's 
3 - Y - 6 method. 

71  - 80:  Integration  method  parameter,  8,  in  Newmark's 

8 - Y - 6 method.  If  zero  or  blank,  8 Is  assumed 
to  be  equal  to  0.25  (1  + 6)2. 
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C2.  STATIC  LOAD  PATTERNS  - NSPAT  sets  of  cards  as  follows. 

Each  set  consists  of  a control  card  followed  by  as  many  cards 
as  needed  to  define  the  nodal  loads.  Load  patterns  are  assumed 
to  be  Input  in  numerical  sequence. 

C2(a)  CONTROL  CARD  (I5,3X,18A4) 

Columns  1 - 5:  Number  of  nodal  load  commands  for  this 

pattern  (NSLC). 

9 - 80:  Load  pattern  title,  to  be  printed  with 
output. 

C2(b)  NODAL  LOADS  (15, 6F10. 0,215)  - NSLC  cards 

Columns  1 - 5:  Node  number,  or  number  of  first  node  in 

a series  of  nodes  covered  by  this  command. 

6 - 15:  Load  In  X-dlrection,  positive  in  positive 
direction  of  X-axis. 

16  - 25:  Load  in  Y-dlrection,  positive  in  positive 
direction  of  Y-axis. 

26  - 35:  Load  in  Z-direction,  positive  in  positive 
direction  of  Z-a\is. 

36  - 45:  Moment  about  X-axis,  positive  by  right  hand  screw 
rule. 

46  - 55:  Moment  about  Y-axis,  positive  by  right  hand 
screw  rule. 

56  - 65:  Moment  about  Z-axIs,  positive  by  right  hand 
screw  rule. 

66  - 70:  Number  of  last  node  in  series.  Leave  blank  for 
a single  node. 

71  - 75:  Node  number  difference  between  successive 
nodes  in  series.  Leave  blank  for  a single 
node,  or  if  node  number  difference  equals 
one. 
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C3.  GROUND  MOTION  (ACCELERATION)  RECORDS 


Omit  if  IGM,  Section  Cl,  Is  zero  or  blank.  Accelerations 
are  assumed  to  be  In  acceleration  units,  not  as  multiples 
of  the  acceleration  due  to  gravity. 

C3(a)  CONTROL  CARD  (4I5.6F10.0)  - One  card 

Columns  1 - 5:  Number  of  time  points  defining  qround  motion 

record  in  X-direction  (NIPX).  Leave  blank 
or  punch  zero  for  no  ground  motion  in  this 
direction. 

6 - 10:  Number  of  time  points  deflnino  ground  motion 
record  in  Y-directlon  (NIPY).  Leave  blank 
or  punch  zero  for  no  ground  motion  in  this 
direction. 


11  - 15: 


Number  of  time  points  defining  qround  motion 
record  in  Z-direction  (NIPZ).  Leave  blank 
or  punch  zero  for  no  ground  motion  in  this 
direction. 


16  - 20: 


21  - 30: 


31  - 40: 

41  - 50: 

51  - 60: 

61  - 70: 

71  - 80: 


Print  code,  as  follows 

a)  zero  or  blank:  records  are  not  printed. 
1:  records  are  printed  as  input  and 
scaled. 

1:  records  are  printed  as  input,  scaled 
and  interpolated  at  time  step  intervals. 


b) 

(c) 


Input  time  interval  for  X-ground  motion.  If 
blank  or  zero,  both  time  and  acceleration  values 
must  be  input;  otherwise  only  acceleration 
values  must  be  input,  the  times  being  auto- 
matically determined.  See  Section  C3(b). 

Input  time  interval  for  Y-qround  motion.  If 
blank  or  zero,  both  time  and  acceleration 
values  must  be  input;  otherwise  only  accelera- 
tion values  must  be  input.  See  Section  C3(c). 

Input  time  interval  for  Z-ground  motion.  If 
blank  or  zero,  both  time  and  acceleration  values 
must  be  input;  otherwise  only  acceleration 
values  must  be  input.  See  Section  C3(d). 

Scale  factor  by  which  X-ground  accelerations 
are  to  be  multiplied. 

Scale  factor  by  which  Y-ground  accelerations 
are  to  be  multipl ied. 

Scale  factor  by  which  Z-ground  accelerations 
are  to  be  multiplied. 
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C3(b)  X RECORD  - One  card  followed  by  as  many  cards  as  needed. 

Omit  If  NIPX  is  blank  or  zero. 

(l)  FIRST  CARD  (15A4.5A4) 

Columns  1 - 60:  Record  title,  to  be  printed  with  output. 

61  - 80:  Input  format  to  read  NIPX  points  defining 
the  record.  For  example,  if  the  format 
is  12F6.0,  punch  (12F6.0). 

(ii)  FOLLOWING  CARDS 

As  many  cards  as  needed  to  specify  NIPX  input  points, 
with  the  format  defined  in  columns  61-80  of  the  first 
card.  If  both  time  and  acceleration  values  are  input, 
the  time  must  imnediately  precede  the  corresponding 
acceleration. 

C3(c)  Y RECORD  - One  card  followed  by  as  many  cards  as  needed. 

Omit  if  NIPY  is  blank  or  zero. 

(1)  FIRST  CARD  (15A4.5A4) 

Columns  1 - 60:  Record  title,  to  be  printed  with  output. 

61  - 80:  Input  format  to  read  NIPY  points  defining 
the  record. 

(ii)  FOLLOWING  CARDS 

As  many  cards  as  needed  to  specify  NIPY  input  points, 
with  the  format  defined  in  columns  61-80  of  the  first 
card. 

C3(d)  Z RECORD  - One  card  followed  by  as  many  cards  as  needed. 

Omit  of  NIPZ  is  blank  or  zero. 

(i)  FIRST  CARD  (15A4.5A4) 

Columns  1 - 60:  Record  title,  to  be  printed  with  output. 

61  - 80:  Input  format  to  read  NIPZ  points  defining 
the  record. 

(ii)  FOLLOWING  CARDS 

As  many  cards  as  needed  to  specify  NIPZ  input  points, 
with  the  format  defined  in  columns  61-80  of  the  first  card. 
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Note:  The  acceleration  scale  factor  may  be  used  to  Increase  or 
decrease  the  accelerations,  or  to  convert  from  multiples 
of  the  acceleration  due  to  gravity  to  acceleration  units. 
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C4. 


C4  (a ) 


C4(b) 


DYNAMIC  FORCE  RECORDS  - NDLGC  sets  of  cards,  as  follows. 

Each  set  consists  of  one  card  followed  by  as  many  cards  as 
needed  to  define  the  record.  Records  are  assumed  to  be  numbered 
in  sequence  as  input. 


FIRST  CARD  (2I5,2Fi0.0,8A4,2X,4A4) 


Columns  1-5:  Number  of  time  points  defining  record 

(NIPT). 

6-10:  Print  code,  as  follows. 

(a)  zero  or  blank:  record  is  not  printed. 

(b)  1:  record  is  printed  as  input  and  scaled. 

(c)  -1:  record  is  printed  as  input  and  scaled 

and  as  interpolated  at  time  step 
intervals. 


11  - 20:  Input  time  Interval.  If  blank  or  zero, 
both  time  and  force  values  must  be  Input; 
otherwise  only  force  values. 

21  - 30:  Scale  factor  by  which  force  values  are  to  be 
mul tipi ied. 

31  - 62:  Record  title,  to  be  printed  with  output. 

65  - 80:  Input  format  to  read  points  defining  the  record. 

FOLLOWING  CARDS 


As  many  cards  as  needed  to  specify  NIPT  Input  points,  with  the 
format  defined  in  columns  65-80  of  the  first  card.  If  both 
time  and  force  values  are  input,  the  time  must  immediately 
precede  the  corresponding  force. 
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C5.  DYNAMIC  FORCE  APPLICATION  (1615)  - NDLGC  Cards  (See  Section  Cl) 


Acceleration  records.  If  specified,  are  applied  automatically, 
assuming  all  support  points  to  move  in  phase.  Force  records 
are  applied  as  defined  by  the  cards  of  this  section. 


Col umns 


1 - 


5:  Dynamic  force  record  number. 

10:  Direction  code,  as  follows. 

(a)  1 

X translation. 

(b)  2 

Y translation. 

(c  3 

Z translation. 

(d)  4 

X rotation. 

(e)  5 

Y rotation. 

(f)  6 

Z rotation. 

Note: 


11  - 80:  Up  to  14  fields,  each  15.  List  the  nodes 

at  which  the  record  is  to  be  applied.  Each 
node  in  the  list  Is  subjected  to  the  scaled 
force  record. 


The  dynamic  forces  as  specified  by  the  dynamic  force  record  number 
are  applied  in  the  positive  direction  defined  by  the  direction  code. 

To  apply  forces  In  the  negative  direction,  the  scale  factor  by 
which  the  force  values  are  multiplied  (Section  C4)  should  be  negative. 


DAMPING  SPECIFICATION  (3F10.0)  - One  card. 

Omit  if  code  for  static  and/or  dynamic  analysis,  KSTAT  (Section  Cl) 
equals  -1. 

Columns  1 - 10:  Mass  proportional  damning  factor,  3^. 

11  - 20:  Tangent  stiffness  proportional  damping 
factor,  See  Note  following. 

21  - 30:  Initial  stiffness  proportional  damping 
factor,  3Q.  See  Note  following. 


If  desired,  It  Is  possible  to  specify  different  values  of  the 
factors  3y  and  3o  for  each  element  group.  See  Section  E for 
explanation  of  this  option. 
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OUTPUT  SPECIFICATION 


This  set  of  cards  consists  of  a control  card  followed  by  as  many 
cards  as  needed  to  specify  node  numbers  for  output.  See  Note 
following. 

a)  CONTROL  CARD  (10I5.7A4)  - One  card. 

Columns  1 - 5 Time  interval  for  printout,  of  nodal  displacement, 

velocity  and  acceleration  time  histories,  expressed 
as  a multiple  of  the  integration  time  step.  Leave 
blank  or  punch  zero  for  no  time  history  output 
or  if  there  is  no  dynamic  analysis. 


6 - 10:  Time  interval  for  printout  of  element  action 

time  histories  (stresses,  forces,  etc.)  expressed 
as  a multiple  of  the  Integration  time  step.  Leave 
blank  or  punch  zero  for  no  time  history  output 
or  If  there  is  no  dynamic  analysis. 

11  - 15:  Time  interval  for  printout  of  Intermediate 
envelopes  of  nodal  displacements  and  element 
actions,  expressed  as  a multiple  of  the  integra- 
tion time  step.  Leave  blank  of  punch  zero  for 
no  intermediate  envelope  output  or  if  there  is 
no  dynamic  analysis.  Envelopes  are  automatically 
output  at  the  end  of  the  dynamic  analysis. 

16  - 20:  Number  of  nodes  for  X-displacement,  velocity 
and  acceleration  output  (NODSC).  For  output 
at  all  nodes,  punch  -1. 

21  - 25:  Number  of  nodes  for  Y-displacement,  velocity 
and  acceleration  output  (NODSY).  For  output 
at  al 1 nodes,  punch  -1 . 

26  - 30:  Number  of  nodes  for  Z-displacement,  velocity 
and  acceleration  output  (NODSZ).  For  output 
at  al 1 nodes,  punch  -1 . 

31  - 35:  Time  interval  for  punched  output  of  nodal 

displacement,  velocity  and  acceleration  time 
histories,  expressed  as  a multiple  of  the 
Integration  time  step.  Leave  blank  or  punch 
zero  for  no  punched  output  or  if  there  is  no 
dynamic  analysis. 

36  - 40:  Number  of  nodes  for  punched  output  of  X- 

displacement,  velocity  and  acceleration  response 
(NODXP) . For  output  at  all  nodes,  punch  -1. 

41  - 45:  Number  of  nodes  for  punched  output  of  Y- 

dlsplacement,  velocity  and  acceleration  response 
(NODYP).  For  output  at  all  nodes,  punch  -1. 


235 


Column  46  - 50:  Number  of  nodes  for  punched  output  of 

Z-displacement,  velocity  and  acceleration 
response  (NODZP).  For  output  at  all  nodes, 
punch  -1 . 

51  - 78:  Format  for  punched  output.  See  Note 
following. 


Note:  Results  for  the  same  nodes  and  elements  are  printed  for  both 

static  and  dynamic  analyses,  except  that  velocities  and  accelera- 
tions are  not  printed  for  static  analyses.  Punched  output  is 
provided  only  for  dynamic  analyses. 

Envelope  values  are  printed  for  the  dynamic  analysis,  and  may  be 
printed  at  the  end  of  each  static  load  increment  if  so  specified 
on  Card  F(a).  For  punched  output,  the  quantities  output  are  node 
number,  direction  ( i . e . , X,  Y or  Z),  displacement,  velocity, 
acceleration  and  time.  The  node  number  and  direction  must  be 
output  in  15  and  A5  format  respectively;  whereas  other  quantities 
may  be  output  in  any  desired  format,  specified  between  parentheses 
in  column  51-78.  For  example 

( I5,A5,3E1 5. 5,1 5X.E10.4) 


D(b)  FOLLOWING  CARDS  - Six  sets  of  cards,  as  follows. 

1.  List  of  nodes  for  X response  printout  (1615)  - As  many  cards 
as  needed  to  specify  NODSX  number  of  nodes,  sixteen  to  a card. 
Omit  if  NODSX  equals  zero  or  -1. 

2.  List  of  nodes  for  Y response  printout  (1615)  - As  many  cards 
as  needed  to  specify  NODSY  number  of  nodes,  sixteen  to  a 
card.  Omit  if  NODSY  equals  zero,  or  -1. 

3.  List  of  nodes  for  Z response  printout  (1615)  - As  many  cards 
as  needed  to  specify  NODSZ  number  of  nodes,  sixteen  to  a card. 
Omit  if  NODSZ  equals  zero,  or  -1. 

4.  List  of  nodes  for  X response  punched  output  (1615)  - As  many 
cards  as  needed  to  specify  NODXP  number  of  nodes,  sixteen  to 
a card.  Omit  if  NODXP  equals  zero,  or  -1. 

5.  List  of  nodes  for  Y response  punched  output  (1615)  - As  many 
cards  as  needed  to  specify  NODYP  number  of  nodes,  sixteen  to 
a card.  Omit  if  NODYP  equals  zero,  or  -1. 

6.  List  of  nodes  for  Z response  punched  output  (1615)  - As  many 
cards  as  needed  to  specify  NODZP  number  of  nodes,  sixteen  to 
a card.  Omit  if  NODZP  equals  zero,  or  -1. 
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E.  ELEMENT  SPECIFICATION 

Elements  must  be  divided  into  "groups. " All  elements  in  any 
groups  must  be  of  the  same  type.  However,  elements  of  the  same 
type  may  be  divided  into  separate  groups  if  desired. 


Elements  groups  may  be  input  in  any 
of  element  groups  may  not  exceed  20. 


sequence.  The  total  number 
The  elements  in  any  group 


must  be  numbered  sequentially,  the  number  of  the  first  element 


in  the  group  being  any  convenient  number. 


El.  THREE  DIMENSIONAL  TRUSS  ELEMENTS 


See  Appendix  R1  for  description  of  element.  Number  of  words 
of  information  per  element  = 96. 


El (a)  CONTROL  INFORMATION  (12I5,4F5.0)  - One  card. 

Columns  5:  Element  group  indicator.  Punch  1 (to  indicate 

that  the  group  consists  of  three  dimensional 
truss  elements). 

6 - 10:  Number  of  elements  in  this  oroun. 

11  - 15:  Element  number  of  the  first  element  in  this 

group.  If  blank  or  zero,  assumed  to  be 
equal  to  1 . 

16  - 20:  Number  of  material  types.  If  blank  or  zero, 
assumed  to  be  equal  to  1 . 

21  - 60:  Blank  (not  used  for  this  element  type). 

61  - 65:  Initial  stiffness  damping  factor  Bo.  If 

blank  or  zero.  Bo  is  assumed  to  be  equal 
to  the  system  Bo  value  input  in  card  C6. 

66  - 70:  Current  tangent  stiffness  damping  factor. 

By.  If  blank  or  zero,  By  is  assumed  to  be 
equal  to  the  system  By  value  input  in  card  C6. 


El (b)  MATERIAL  PROPERTY  INFORMATION  (I5.4F10.0)  - One  card  for  each 
different  material  type. 

Columns  1 - 5:  Material  number,  in  sequence  starting  with  1. 

6-15:  Young's  modulus  of  elasticity,  E. 

16  - 25:  Strain  hardening  modulus  as  a proportion 
of  Youno's  modulus  (i.e.,  the  ratio  E^/E) . 
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Column  26  - 35: 


Yield  stress  in  tension. 


36  - 45:  Yield  stress  In  compression,  or  elastic 

buckling  stress  in  compression  (Input  as 
a positive  value). 


ELEMENT  GENERATION  COMMANDS  (415, 2F10. 0,415)  - As  many 
cards  as  needed  to  generate  all  elements  in  this  group. 


Cards  must  be  entered  in  order  of  increasing  element  number. 
Cards  for  the  first  and  last  element  must  be  included.  See 
Note  for  explanation  of  generation  procedure. 


Columns  1 - 5:  Element  number,  or  number  of  first  element 

in  a sequentially  numbered  series  of  elements 
to  be  generated  by  this  card. 


6 - 10:  Node  number  at  element  end  i. 

11  - 15:  Node  number  at  element  end  j. 

16  - 20:  Material  number.  If  blank  or  zero,  assumed 

to  be  equal  to  1 . 

21  - 30:  Cross-sectional  area. 

31  - 40:  Initial  axial  force  on  the  element. 

41  - 45:  Node  number  increment  for  element  generation. 

If  blank  or  zero  assumed  to  be  equal  to  1. 

50:  Code  for  large  displacement  effects.  Leave 

blank  or  punch  zero,  for  small  displacement 
effects.  Punch  1 for  large  displacement 
effects. 

55:  Time  history  output  code.  Leave  blank  or 

punch  zero  for  no  time  history  output. 

Punch  1 if  time  history  output  is  required. 

60:  Buckling  code.  Leave  blank  or  punch  zero  if 

element  yields  in  compression  without  buckling. 
Punch  1 if  element  buckles  elastically  In 
compression. 
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E4.  THREE  DIMENSIONAL  EIGHT  TO  TWENTY  NODE  SOLIDS 
E4(a)  CONTROL  INFORMATION  (12I5,4F5.0)  - One  card. 


> 


Columns 


1 - 5:  Element  group  indicator.  Punch  4 (to  indicate 

that  the  grouo  consists  of  three-dimensional 
solid  elements). 

6 - 10:  Number  of  elements  In  this  group. 

11  - 15:  Element  number  of  the  first  element  in 

this  group. 

16  - 20:  Number  of  nodes  (8  to  20  nodes  per  element 
can  be  specified). 

21  - 25:  Number  of  material  types. 

26  - 30:  Material  model  number  as  follows: 


(a)  1:  Linear  elastic  isotropic  model. 

(b)  2:  Linear  elastic  orthotropic  model. 

(c)  3:  Elastic  perfectly  nlastlc.  Von 

Misces  yield  criteria. 

(d)  4:  Ramberg-Osgood  model. 

(e)  5:  Elastic  Plastic  Fracture  strain- 

hardening  model. 


31  - 

35: 

r-direction  integration  order. 

36  - 

40: 

s-direction  Integration  order. 

41  - 

45: 

t-direction  Integration  order. 

61  - 

65: 

Initial  stiffness  damping  factor. 

66  - 

70: 

Tangent  stiffness  damping  factor. 

E4(b)  MATERIAL  INFORMATION 


E4(b) . 5 ELASTIC  PLASTIC  FRACTURE  STRAIN  HARDENING  MATERIAL  MODEL 


E4(b) . 5. 1 GENERAL  MATERIAL  INFORMATION  (I5,2F10.0,F5.0,2F10.0)  - One  card. 


Columns 

1 - 

5: 

Material  Identification  number, 
indicate  that  model  is  elastic  t 
strain  hardening  material  model) 

Punch  5 (to 
ilastlc  fracture 

i . 

6 - 

15: 

Young's  Modulus  of  Elasticity  E. 

16  - 

25: 

Poisson's  ratio. 

26  - 

30: 

Number  of  Increments  describing  the  nonelastic 
portion  of  uniaxial  stress-strain  curve  (KKCC). 
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31  - 40:  Coefficient  by  which  the  Ce  matrix  is  to  be 
multiplied  after  cracking  or  crushing  of 
concrete.  Punch  1 In  column  40  if  material 
is  not  concrete. 

E4(b) . 5. 2 ELASTIC  PLASTIC  FRACTURE  STRAIN  HARDENING  MATERIAL  INFORMATION 
(4F5.0.6F10.0)  - One  card. 

Columns  1 - 5:  K'-  value  for  compression  region  (default 

value  = 3). 

6 - 10:  KJ-value  for  tension  region  (default 
value  = 3). 

11  - 15:  X-value  for  compression  region  (default 
value  = 3) . 

16  - 20:  X-value  for  tension  reqion  (default  value 

= -3). 

21  - 30:  Ultimate  compressive  stress  (must  be 
greater  than  zero) . 

31  - 40:  Ultimate  tensile  stress/ul timate  compressive 
stress  (must  be  greater  than  zero). 

41  - 50:  Ultimate  biaxial  comnresslve  stress/ultimate 
compressive  stress. 

51  - 60:  Yield  comoressive  stress/ultimate  compressive 
stress. 
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E4(c)  ELEMENT  GENERATION  COMMANDS  - Two  cards  per  element  If  the  number 
of  nodes  per  element  is  greater  than  8.  As  many  sets  of  these 
two  cards  as  needed  to  generate  all  the  elements  in  this  group. 

E4(c) . 1 ELEMENT  CARD  (13I5.3F5.0) 

Columns  1 - 5:  Element  number. 

6 - 10:  Node  number  at  element  node  1. 

11  - 15:  Node  number  at  element  node  2. 

16  - 20:  Node  number  at  element  node  3. 

21  - 25:  Node  number  at  element  node  4. 

26  - 30:  Node  number  at  element  node  5. 

31  - 35:  Node  number  at  element  node  6. 

36  - 40:  Node  number  at  element  node  7. 

41  - 45:  Node  number  at  element  node  8. 

46  - 50:  Material  model  number. 

51  - 55:  Node  number  increment  for  element 
generation. 

56  - 60:  Large  displacement  effects  code. 

(a)  1:  include  large  displacement 

effects. 

(b)  0:  small  displacement  effects 

only. 

61  - 65:  Time  history  output  code,  printout 
of  stresses  and  strains  as  follows: 

(a)  0:  no  printout. 

(b)  1:  printout  at  output  point. 

(c)  2:  printout  at  gauss  points 

and  output  point. 

66  - 70:  r-coordinate  of  output  point. 

71  - 75:  s-coordinate  of  output  point. 

76  - 80:  t-coordinate  of  output  point. 
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I 

E4(c) . 2 ADDITIONAL  NODES  SPECIFICATION  (1215)  - Omit  this  card  if 
number  of  nodes  is  less  than  or  equal  to  8. 

Columns  1 - 5:  Node  number  at  element  node  9. 

6 - 10:  Node  number  at  element  node  10. 

II  - 15:  Node  number  at  element  node  11. 

16  - 20:  Node  number  at  element  node  12. 

21  - 25:  Node  number  at  element  node  13. 

26  - 30:  Node  number  at  element  node  14. 

31  - 35:  Node  number  at  element  node  15. 

36  - 40:  Node  number  at  element  node  16. 

41  - 45:  Node  number  at  element  node  17. 

46  - 50:  Node  number  at  element  node  18. 

51  - 55:  Node  number  at  element  node  19. 

56  - 60:  Node  number  at  element  node  20. 
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F.  STATIC  ANALYSIS  SPECIFICATION  - NSLGC  sets  of  cards  (see  Section  Cl). 

Each  set  consists  of  a solution  procedure  card  followed  by  one 
or  more  cards  defining  a linear  combination  of  static  force  patterns. 
Each  set  defines  an  increment  of  static  load. 


F(s). SOLUTION  PROCEDURE  CARD  (815,  4F10.0)  - One  Card 

Columns  1 - 5:  Number  of  equal  steps  in  which  load  increment 

is  to  be  applied,  positive  if  results  envelopes 
are  not  to  be  printed  at  the  end  of  the  increment, 
otherwise  negative. 

6 - 10:  Iteration  type,  as  follows: 

(a)  zero  or  blank:  Newton-Raphson  iteration 

(b)  n:  Constant  stiffness  iteration  with 

alpha-constant  over-relaxation,  the 
alpha  matrix  being  reinitialized  every  n 
iterations. 

15:  Type  of  state  determination  calculation  to  be 
used  for  constant  stiffness  iteration  as 
f ol 1 ows : 

(a)  zero  or  blank:  path  independent. 

(b)  1:  path  dependent. 

Path  dependent  state  determination  is  always 
used  for  Newton-Raphson  iteration. 

16  - 20:  Stiffness  reformation  code,  as  follows. 

(a)  zero  or  blank:  stiffness  used  in 

preceding  step  is 
retained. 

(b)  n:  stiffness  is  reformed  every  n 

load  steps. 

25:  Termination  code,  as  follows. 

(a)  zero  or  blank:  If  the  solution  does  not 

converge  within  the  maximum 
number  of  iterations  for 
any  load  step,  the  next 
load  step  will  be  applied. 

(b)  1:  If  the  solution  does  not  converge, 

the  execution  will  terminate. 

26  - 30:  Print  code,  as  follows. 

(a)  -1:  results  are  not  printed  for  this 

increment. 

(b)  zero  or  blank:  results  are  printed  at 

the  end  of  the  increment 
only. 

(c)  1:  results  are  printed  after  each  load  step. 

(d)  2:  results  are  printed  every  iteration. 

This  option  should  be  used  for  debugging 
purposes  only. 
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Columns  31  - 35:  Maximum  number  of  cycles  of  iteration  within 

any  load  step.  If  set  to  -1,  the  unbalanced 
load  will  be  set  to  zero  at  a new  load  step- 
works  only  for  imposed  displacements  option. 

36  - 40:  Maximum  number  of  iterations  within  any  cycle. 

41  - 50:  Nodal  force  convergence  tolerance  to  be  used 

in  last  step  of  load  increment. 

51  - 60:  Nodal  force  convergence  tolerance  to  be  used 
in  all  except  last  step  of  load  increment. 

61  - 70:  Nodal  force  tolerance  for  change  of  stiffness 
in  Newton-Raphson  iteration.  If  the  unbalanced 
force  reduces  below  this  tolerance,  the  stiff- 
ness will  not  be  reformed  for  the  next 
iteration. 

71  - 80:  Maximum  nodal  displacement  (translation  or 

rotation)  increment  permitted  in  any  iteration 
step.  Leave  blank  for  unlimited  displacement. 
Displacement  limits  should  be  specified  only 
with  Newton-Raphson  iteration. 

F(b).  FOLLOWING  CARDS  (8F10.0)  - As  many  cards  as  needed. 

Columns  1 - 80:  Up  to  eight  fields,  each  F10.0.  For  each 

static  force  pattern  in  turn,  specify  a scale 
factor  by  which  the  Dattern  is  to  be  multi- 
plied. The  scaled  patterns  are  added  together 
to  produce  the  load  increment. 

Scale  factors  may  be  positive  or  negative. 

Leave  the  corresponding  field  blank  or 
punch  zero  to  ignore  anv  force  nattern. 

F(c)  IMPOSED  DISPLACEMENT  CONTROL  (I5.F10.0)  - Omit  this  card  if  NC0RD  = 0 
(See  section  B1 ) - One  card. 

Columns  1 - 5:  Imposed  displacements  identification  as 

follows : 

(a)  0:  no  displacements  imposed. 

(b)  1:  imposed  displacements. 

6-15:  Scale  factor  by  which  the  imposed  displacement 

pattern  is  to  be  multiplied.  The  scale  factor 
multiplied  to  the  imposed  displacement  pattern 
is  the  net  value  of  imposed  displacements  in  a 
single  load  step. 
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NEW  PROBLEM 


Data  for  a new  problem  may  follow  1 mned lately  startlnq  with 
Section  A.  Any  number  of  structures  may  be  analyzed  In  a 
single  computer  run. 


I.*_ TERMINATION  CARD  (A4)  - One  card  to  terminate  the  complete  data  deck. 

Columns  1 - 4:  Punch  the  word  STOP. 


i 


245 


2 


SAMPLE  INPUT  TO  PROGRAM  FOR  CCS-?  (REFER  FIGURE  45) 


START  * CCS-2  * RHO  ■ 2 PER  CENT  * PURE  COMPRESSION  * 
40  20  0 2 0 0 2 0 0 -1 


1 

314.96 

4 

354.33 

5 

222.71 

222.71 

8 

250.55 

250.55 

9 

314.96 

12 

354.33 

13 

. 

222.71 

222.71 

16 

- 

250.55 

250.55 

17 

■ 

314.96 

20 

- 

354.33 

21 

39.37 

314.96 

24 

39.37 

354.33 

25 

39.37 

222.71 

222.71 

28 

39.37 

250.55 

250.55 

29 

39.37 

314.96 

32 

39.37 

354.33 

33 

39.37 

222.71 

222.71 

36 

39.37 

250.55 

250.55 

37 

39.37 

314.96 

40 

39.37 

354.33 

1 

000111 

40 

1 

111111 

20 

1 

2 1 

4 -1.0 

0.0 

29 

32  1 

21 

24  1 

25 

28  1 

33 

36  1 

37 

40  1 

1 

1 2 

1 

1 

CCS-2 

PURE  COMPRESSION 

-1  -1 

-1 

1 

26 

1 

1 

1 

30000. 

0.02 

60.0 

60.0 

1 

1 

21 

1 

82.78 

3 

0 

1 

0 

3 

5 

25 

1 

82.78 

3 

0 

1 

0 

5 

9 

29 

1 

82.78 

3 

0 

1 

0 

7 

13 

33 

1 

82.78 

3 

0 

J 

0 

9 

17 

37 

1 

82.78 

3 

0 

1 

0 

11 

1 

5 

1 

15.50 

4 

0 

1 

0 

15 

4 

8 

1 

15.50 

4 

0 

1 

0 

19 

21 

25 

1 

15.50 

4 

0 

1 

0 

23 

24 

28 

1 

15.50 

4 

0 

1 

0 

26 

36 

40 

1 

15.50 

4 

0 

1 

0 

4 

12 

1 

8 

1 

5 

2 

2 

2 

1 

4000. 

000000001 

4 

0.10 

4.0 

5.0 

0.12 

1.15 

0.48 

0.10 

0.552 

0.48 

0.0006 

0.717 

0.001 

1.0 

0.002 

0.65 

0.003 

1 

4 

8 

28 

24 

3 

7 

27 

23 

1 

-1 

2 

4 

8 

12 

32 

28 

7 

11 

31 

27 

1 

-1 

0 

2 

7 

13 

9 

29 

33 

14 

10 

30 

34 

1 

1 

0 

2 

10 

17 

13 

33 

37  1 

18 

14 

34 

38 

1 

1 

0 

2 

12 

19 

15 

35 

39  20 

16 

36 

40 

1 

1 

0 

2 

2 

1 

0 

1 

0 

1 

5 

2 

1.0 

1.0 

1.0 

1 

0.02 

5 

1 

0 

1 

0 

1 

5 

2 

1.0 

1.0 

1 

.0 
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APPENDIX  B 


1.  FIBER  ANALYSIS:  BALANCED  POINT  ON  P-M  INTERACTION  CURVE 

Assume  the  strain  distribution  and  stress  distribution  as  shown  In 
Figure  B.l.  The  stress-strain  curves  for  concrete  and  steel  are  shown 
in  Figures  17  and  18,  respectively. 

a.  Forces  in  the  Tension  and  Compression  Blocks 
Since  steel  undergoes  yield 

Ts  ■ Vs  0) 

T ■ 60  x 31  * 1860  kips 


Compressive  steel  force: 


C$  = 30000  x 0.00129  x 31 
C$  = 1203  kips 


Compressive  block  of  concrete: 


m 


c - r be  [ - 

c c L c 


co 


-1( 

3 1 cco 


)2  ] 


(2) 


(3) 


where 

b = width  of  the  section 

f*  = ultimate  stress  of  concrete  in  uniaxial  compression 
c 

c = depth  of  the  neutral  axis  from  the  top  fiber 
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Interaction  Curve 


e?  = strain  in  extreme  fiber  of  concrete 
e = strain  in  concrete  at  maximum  stress 

Cc  = 5.0  x 39.37  x 16.7325  [1  - 1/3] 

Cc  = 2195.86  kips 


b.  Calculation  of  P and  M: 


Cc  + Cs  = P + Ts  (4) 

P = 2195.86  + (203  - 1860 
P = 1538  kips 

Moments  about  centroid  of  section: 

Centroid  of  compress  ion  center  y is  given  by 


y = 


c[2/3  (c;/Eco)  - 1/4  (e^co)-] 
- '/3  <#««.>* 


(5) 


y=  ,6.7325  L2/3rr_^41 

y = 10.46  ins. 

M = Ccy  + Csy;  + T$ys  (6) 

where 

y = centroid  of  concrete  compressive  block 
y^  = lever  arm  to  the  centroid  of  section  of  compression 
steel 
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r 


y$  = lever  arm  to  the  centroid  of  section  of  tension  steel 
M = 2195.86  x 13.4125  + 1538  x 13.78  + 1860  x 13.78 
M = 76.276.41  kip-ins 

c.  Calculation  of  deformations  6 and  0 at  the  centroid  of  the  section: 


f cent  = 0.000353  in/in 
^ cent  ~ ^ent*- 
where 

15 cent  = disPlacement  at  centroid 
= strain  at  centroid 

cent 

L = length  of  section 

5 cent  = 0-000353  x 19.685 

S„.  = 0.00695  ins 
cent 

0 4.  = 4.  L 

cent  Ncent 

where 

0cent  = rotation  at  geometric  centroid 
<|>cent  = curvature  at  geometric  centroid 

4 = 0.002  - 0.000353 

9cent  13.78 

-0Cent  = 0.0001195  l/1n 

0Cent  = 0.0001195  x 19.685 

0 . = 0.002353  rad. 

cent 


(7) 


(8) 


» 
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d.  Ratio  of  5 to  0 


6/0  = 0.00695/0.002353 
5/0  = 2.95 

Note  that  the  displacement  6 Is  an  extension  displacement 
at  the  centroid. 
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Los  Alamos  Scientific  Laboratory 
P.0.  Box  1663,  Los  Alamos,  NM  87545 
1 ATTN:  Doc  Control  for  Reports  Lib 

Sandia  Laboratories 

Livermore  laboratory,  P.0.  Box  969,  Livermore,  CA 
1 ATTN:  Doc  Control  for  Tech  Lib 

Sandia  Laboratories 

P.0.  Box  5800,  Albuquerque,  NM  87115 
1 ATTN:  Doc  Con  for  3141  Sandia  Rpt  Coll 

U.S.  Ehergy  Rsch  & Dev  Admin,  Div  of  Hq  Services 
Library  Branch  G-043 
Washington,  DC  20545 

1 ATTN:  Doc  Con  for  class  Tech  Lib 

U.S.  Ehergy  Rsch  & Dev  Admin 

Nevada  Opeaations  Office,  P.0.  Box  5400,  Albuq,  NM 
1 ATTN:  Doc  for  Tech  Lib 

U.S.  Ehergy  Rsch  & Dev  Admin,  Div  of  Hq  Services 
Library  Branch  G-043 
Washington,  DC  20545 

1 ATTN:  Doc  Con  for  Class  Tech  Lib 

US  Ehergy  Rsch  & Dev  Admin 
Nevada  Operations  Office 
P.0.  Box  14100,  Las  Vegas,  NV  89114 
1 ATTN:  Doc  Con  for  Tech  Lib 

Dept  of  the  Interior 
Bureau  of  Mines 

Bldg  20,  Denver  Federal  Center,  Denver,  CO  80225 
1 ATTN:  Tech  Lib 

Department  of  Defense  Contractors 

Aerospace  Corporation 
P.0.  Box  92957,  Los  Angeles,  CA  90009 
1 ATTN:  Tech  Info  Services 

Agbabian  Associates 

250  N.  Nash  St.,  El  Segundo , CA  90245 
1 ATTN:  M.  Agbabian 

Analytic  Services,  Inc. 

5613  Leesburg  Pike,  Falls  Church,  VA  22041 
1 ATTN:  George  Hesselbacher 


94550 


87115 
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Applied  Theory,  Inc. 

1010  Westwood  Blvd , Los  Angeles,  CA  90024 
1 ATTN:  Jolm  C.  Trullo 

Artec  Associates,  Inc. 

2 (>046  Eden  land  ing  Road,  I Lay  ward  CA  94545 
1 ATTN:  Steven  Gill 

Avco  Research  & Systems  Group 

201  Lowell  Street,  Wilmington  HA  0^887 
1 ATTN:  Research  Lib  A«30  Rm  7201 

lkittelle  Memorial  Institute 

505  King  Ave. , Columbus  OH  43201 
1 ATTN:  Tech  Lib 

RDM  Corporat  Lon 

1920  Aline  Ave,.  Vienna  VA  22180 
1 ATTN:  Tech  Lib 

Ted  Rely t schko 

b304  No.  Hiawatha  Ave. , Vienna  VA  22180 
1 ATTN:  Ted  Bely t schko 

Boeing  Company 

P . 0 . Box  3607,  Seattle,  WA  98124 
1 ATTN:  Aerospace  Lib 

Brown  Engineering  Co. , Inc 

Cummings  Research  Park,  Huntsville  AL  35807 
L ATTN:  Manu  Patel 

California  Institute  of  Tec  lino  logy 

1201  E.  California  Blvd,  Pasadena  CA  91109 
1 ATTN:  Thomas  J.  Alvrens 

California  Research  & Technology,  Inc. 

6269  Variel  Ave.,  Woodland  Hills  CA  91364 
1 ATTN:  Tech  Lib 

Calspan  Corporation 

P.0.  Box  235,  Buffalo  NY  14221 
1 ATTN:  Tech  Lib 
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Civil/Nuclear  Systems  Corp. 

1200  University  Blvd. , NE. , Albuquerque  NM  87102 
ATTN:  Tech  Lib 

Dayton,  University  of 
Industrial  Security  Super  KL-505 
300  College  Park  Ave. , Dayton,  OH  45409 
ATTN:  Hal lock  F.  Swift 

University  of  Denver 

Colorado  Seminary 

Denver  Research  Institute 

p.O.  Box  10127,  Denver,  CO  80210 

ATTN:  Sec  Officer  for  Tech  Lib 

HG&G,  Inc. 

Albuquerque.  Division 

P.O.  Box  10218,  Albuquerque,  NM  87114 

ATTN:  Tech  Lib 

Electric  Power  Research  Institute 
3412  Hillview  Ave.,  Falo  Alto,  CA  94303 
ATTN:  George  Sifter 

Engrg  Decision  Analysis  Co.  , Inc. 

2400  Michelson  Dr.,  Irvine,  CA  92715 
ATTN:  R.  P.  Kennedy 

Franklin  Institute 

20th  Street  and  Parkway,  Philadelphia,  PA  19103 
ATTN:  Zen  on  s Zudnns 

General  Electric  Co.,  TIHPO-Centor  for  Advanced  Studies 
8 1 (>  State  St.,  (P.O.  Drawer  i^5)»  Santa  Barbara,  CA  93102 
ATTN:  DASIAC 

General  Research  Corporation 

P.O.  Box  3587,  Santa  Barbara,  CA  93105 

ATTN:  Benjamin  Alexander 

H-Tech  Laboratories,  Ind. 

P.O.  Box  lb8(>,  Santa  Monica,  CA  90406 
ATTN:  B.  Hartenbaum 
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IIT  Research  Institute 
10  West  35th  St.,  Chicago,  IL  60616 
1 ATTN:  Tech  Lib 


Institute  for  Defense  Analyses 

400  Army-Navy  Dr.,  Arlington,  VA  22202 
1 ATTN:  IDA  Librarian  Ruth  S.  Smith 

J.  H.  Wiggins,  co. , Inc. 

1650  S.  Pacific  Coast  Highway,  Redondo  Beach,  CA  90277 
1 ATTN:  John  Collins 

Kaman  Avidyne,  Division  of  Kaman  Sciences  Corp. 

83  Second  Ave. , Northwest  Industrial  park,  Burlington,  MA 
1 ATTN:  Tech  Lib 

Kaman  Sciences  Corporation 

P.0.  Box  7463,  Colorado  Springs,  CO  80933 
l ATTN:  Library 

Karagozian  and  Case 

6330  N.  Figueroa  St.,  Los  Angeles,  CA  90042 
1 ATTN:  Jolin  Karagozian 

Lockheed  Missiles  & Space  Co.  , Inc. 

P.0.  Box  504,  Sunnyvale,  CA  94088 
1 ATTN:  Tech  Lib 

Martin  Marietta  Aerospace,  Orlando  Div 

P.0.  Box  5836,  Orlando,  FL  32805 
1 ATTN:  G.  Fotieo 

McDonnell  Douglas  Corp. 

5301  Bo Isa  Ave.,  Huntington  Beach,  CA  92647 
1 ATTN:  Robert  W.  Halprin 

Merritt  Cases,  Inc. 

*.0.  Box  1206,  Redlands,  CA  92373 
1 ATTN:  Tech  Lib 

Newmark,  Nathan  M. , Consulting  Engineering  Services 

1211  Civil  Engineering  Bldg.,  Rm  B106A,  Hniv  of  Illinois 
1 ATTN:  Nathan  M.  Newmark 
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TRW  Systems  Group 

One  Space  Park,  Redondo  Beach,  CA  90268 
1 ATTN:  Tech  Info  Center/S-1930 

1 ATTN:  Peter  K.  Dai,  Rl/2170 

1 ATTN:  Norm  Lipner 

TRW  Systems  Group,  San  Bernardino  Operations 
P.0.  Box  1310,  San  Bernardino,  CA  92402 
1 ATTN:  E.  Y.  Wong,  527/712 

TRW  Systems  Group 

Room  712,  Bldg  572, Norton  AFB,  CA  92409 
1 ATTN : Gregory  D.  Hulcher 

Universal  Analytics,  Inc. 

7740  W.  Manchester  Bid. , Playa  Del  Rey,  CA  90291 
1 ATTN:  E.  I.  Field 

URS  Research  Co. 

155  Bovet  Rd.,  San  Mateo,  CA  94402 

1 ATTN:  Tech  Lib 

Eric  H.  Wang  Civil  Ehgineering  Rsch  Fa,  University  Station 
Box  188,  University  of  NM,  Albuquerque,  NM  87121 

2 ATTN:  Jerry  Berglund 

Weidlinger  Assoc.  Consulting  Ehgineers 
110  East  59th  St.,  New  York,  NY  10022 
1 ATTN:  Melvin  L.  Baron 

Weidlinger  Associates  Consulting  Engineers 
Suite  245,  3000  Sand  Hill  Road,  Menlo  Park,  CA  94025 
1 ATTN:  F.  S.  Wong 

1 ATTN:  J.  Isenberg 

Dept  of  Civil  Engineering 
University  of  California 
1 ATTN:  Prof  Alex.  C.  Scordelis 

729  Davis  Hall 
Berkeley,  CA  94720 

Purdue  University 
School  of  Civil  Engineering 
West  Lafayette,  IN  47907 
1 ATTN:  Prof  W.  F.  Chen 
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1 


1 


1 


Pacifica  Technology 

P.0.  Box  148,  Del  Mar,  CA  92014 

ATTN:  G.  Kent 

ATTN:  Prof.  Robert  Dunham 

Physics  International  Co. 

2700  Merced  St. , San  Leandro,  CA  94577 
ATTN:  Doc  Con  for  Tech  Lib 

R & D Associates 

P.0.  Box  9695,  Marine  Del  Rey,  CA  90291 
ATTN:  Tech  Lib 

Rand  Corporation 

1700  Main  St.,  Santa  Monica,  CA  90406 
ATTN:  Tech  Lib 

Science  Applications,  Inc. 

2201  San  Pedro  NE,  Albuquerque,  NM  87110 
ATTN:  J.  L.  Bratton 

Science  Applications,  Inc. 

P.  0.  Box  2351,  La  Jolla,  CA  92038 
ATTN:  Tech  Lib 

Southwest  Research  Institute 

P.0.  Box  Drawer  28510,  San  Antonio,  TX  78284 

ATTN:  Wilfred  E.  Baker 

Stanford  Research  Institute 

333  Ravenswood  Ave. , Menlo  Park,  CA  94026 

ATTN:  Carl  Peterson 

Systems,  Science  and  Software,  Inc. 
p.O.  Box  1620,  La  Jolla,  CA  92038 
ATTN:  Tech  Lib 
ATTN:  G.  Hegemeir,  Room  2091 

Terra  Tek,  Inc. 

420  Wakara  Way,  Salt  Lake  City,  UT  87108 
ATTN:  Tech  Lib 

Tetra  Tech,  Inc 

630  North  Rosemead  Blvd. , Pasadena,  CA  91107 
ATTN:  Tech  Lib 
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PMB  Systems  Ehg 
500  Sansome  Street 
San  Francisco  94111 
ATTN:  Dr.  R.  Litton 


15 


